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ABSTRACT 

We presented the results of an analysis of four XMM-Newton observations of the 
starburst galaxy IC342 taken over a four-year span from 2001 to 2005, with an em- 
phasis on investigating the long-term flux and spectral variability of the X-ray point 
sources. We detected a total of 61 X-ray sources within 35' x 30' of the galaxy down 
to a luminosity of (1-2) x 10 37 erg s -1 depending on the local background. We found 
that 39 of the 61 detected sources showed long-term variability, in which 26 of them 
were classified as X-ray transients. We also found 19 sources exhibiting variations in 
hardness ratios or undergoing spectral transitions among observations, and were iden- 
tified as spectral variables. In particular, 8 of the identified X-ray transients showed 
spectral variability in addition to flux variability. The diverse patterns of variability 
observed is indicative of a population of X-ray binaries. We used X-ray colors, flux 
and spectral variability, and in some cases the optical or radio counterparts to classify 
the detected X-ray sources into several stellar populations. We identified a total of 11 
foreground stars, 1 supersoft sources (SSS), 3 quasisoft sources (QSS), and 2 super- 
nova remnants (SNR). The identified SSS/QSS are located near or on the spiral arms, 
associate with young stellar populations; the 2 SNR are very close to the starburst 
nucleus where current star formation activities are dominated. We also discovered a 
spectral change in the nuclear source of IC342 for the first time by a series of X-ray 
spectrum analysis. 

Subject headings: galaxies: individual (IC342) — X-rays: galaxies — galaxies: star- 
burst 



1. INTRODUCTION 



IC342 is a nearby (1.8Mpc; see iButa & McCalllll999l for a review) late-type Sc/Scd galaxy 
in the Maffei Group which is one of the close st groups to ou r Galaxy. Its spiral arms are well 



developed and are almost face-on (i = 25° ± 3° : iNewtonll 19801) . Its active star formation activities 



in th e nuclear region has made it a popu l ar target for infrared a nd submillimeter observations 
(e.g. iBoker. van der Marel. & Vaccal Il999l Meier & Turner! 120051) . These observations revealed 
that the physical properties of molecular clouds, the infrared luminosity, and the presence of a 



nuclear stellar clus ter in IC342 are similar to the Milky Way in many ways (|Schinnerer et al.ll2003 



Schulz et al.ll200l|) . Its proximity and its orientation provide a unique possibility to study the X- 
ray sources in IC342 from a very favorable viewing angle. One major drawback of X-ray study 
is that this galaxy is located near the galactic plane (b =10.6° ) and has a large absorption (Nn = 
3 x 10 21 cm -2 ) towards the Galactic center ( Stark et al. 1992|) . This limits us to constrain local 
absorption and determine X-ray emission below 1 keV. On the other hand, IC342 has been studied 
in the X-ray with many missions since Einstein. With the advent of high angular resolution and 
high sensitivity instruments such as Chandra and XMM-Newton it is possible to study the X-ray 
source population in depth. 



The first X-ray observation of IC342 from Einstein (|Fabbiano & Trinchierilll987l) showed 
that the X-ray emission was dominated by three ultraluminous X-ray sources (ULX s) which had 



lumin osities above 10 39 erg s 1 (named X-l, X-2, and X-3 based on designations of lOkada et al. 



19981) . Subsequent ROSAT HRI observation in 1991 revealed seven additional poi nt sources in 
the disk with a detection limit of L0.1-2.5keV ~ 2 x 10 37 er gs" 1 dB regma n et alJI 1993b. Two ASCA 



observations taken in 1993 and 2000 on X-l and X-2 (Okadaet al 



1998 



Kubota et all 2002) 



showed that both X-l and X-2 exhibited spectral and inten sity t r ansitio ns w hich resemb l ed the 
spectral/intensity states of X-ray transients in our Galaxy. llCong ( 2003 ) and Bauer et al. ( 2003 ) 
presented the result from a 2001 XMM-Newton 10 ks observation and both reported detections 
of about 35 sources. The slope of the X-ray luminosity function was found to be ~ 0.5 in both 
studies, suggesting the X-ray population of IC342 was consistent with other starburst galaxies 
and Galactic HMXBs. A de tailed analysis on the Chandra HRC-I image observed in 2006 was 
presented in iMak et al] (120081) in which the nuclear X-ray source (X-3) was found to be spatially 
coincident with a nuclear stellar clu ster which had been identified in earlier optical and infrared 
observations (e.g. IBoker et alJll999l) . Thus they proposed that the source was not an ULX and 
was instead associated with starburst activities, together with the possibility of an embedded low 
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luminosity AGN in the nuclear stellar cluster. 



X-ray observations of our own Galaxy revealed a diverse population of sources, including low- 
mass X-ray binaries (LMXBs), high-mass X-ray binarie s (HMXBs), and a few young supernova 



remnants (SNRs) that show various kinds of variability (Fa bbianoll2006|) . Variability on flux and 
spectral properties are useful for characterization of the emission mechanism of X-ray sources. 
However, not much work had been done on the long-term variability of the X-ray sources in IC342 
though there were multiple epochs over the previous two decades. Recently IC342 was observed 
four times with XMM-Newton between 2001 and 2005, and five times with Chandra between 2002 
and 2006. Therefore, we made use of archival XMM-Newton observations which spanned a period 
of four years to study the nature of the X-ray population in IC342. 

In this paper, we present the source catalog and source properties in the four XMM-Newton 
observations. The observation and the data reduction procedures of the X-ray data are presented 
in Section [2] In Section [3j the analysis of the data, including source detection, photometry, com- 
putation of hardness ratios of the X-ray sources are described. In Section HI we present the results 
of the identification and variability properties of the detected X-ray sources. We finally discuss the 
global X-ray source population in IC342, with detailed descriptions of several bright individual 
sources in Section \5\ 



2. OBSERVATIONS AND DATA REDUCTION 

IC342 was observed with the European Photon Imaging Camera (EPIC) on board XMM- 
Newton four times from 2001 to 2005. The instrument modes were full-frame with medium optical 
blocking filter for the EPIC-PN and the two EPIC-MOS camera. The observation in 2001 was part 
of a multi-wavelength campaign to study star formation activities of IC 342 and the aimpoint was at 
the nuclear region. The three subsequent observations were made to study the ULXs in IC 342 with 
two long observations in 2004 (r exp > 20 ks) aiming at X-l and one short observation (r exp w 6 ks) 
aiming at X-2 in 2005. These observations covered the entire galaxy with a circular field of view 
of 30' diameter. Since the field of view of these observations overlapped and were made at various 
aimpoints, the total region covered was about 35' x 30'. Therefore sources at the outer edges of the 
CCD might not be observable in all four exposures. A summary of these observations was listed 



in Ta ble [p E ven though the 2 001 February observation had already been studied in detail by lKong 



J2003I) and lBauer et al.l (|2003h . for for consistency we included it also in our present analysis and 



reanalysed it whenever necessary. 

The event files were filtered and reprocessed using the XMM-Newton Sciences Analysis Sys- 
tem (SAS v7.1.0). Time intervals contaminated by soft-proton flares were identified using the 
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background light curves in the > 10 keV band. High background level with count rate over 0.2 cps 
for MOS and 1 .0 cps for EPIC-PN were excluded and the good time intervals (GTI) were obtained 
for each observation which ranged from 5.6 to 23.6 ks. The resulting GTI of the two observations 
in 2004, with exposure time over > 20 ks for each MOS camera, roughly double those of the 2001 
and 2005 data. Therefore the sensitivity of the 2004 observations were expected to be higher and 
fainter sources could be detected. We selected only good event patterns for imaging: < 12 for 
MOS and < 4 for PN, and restricted our analysis in the energy range 0.2-12 keV. These filtered 
event files were then used for data analysis using HEAsoft v6.4 and XSPEC vl2.4. 



3. ANALYSIS 



3.1. Source Detection and Astrometry 



Source detection using maximum likelihood fitting was done simultaneously on each of the 
EPIC-PN, MOS1, and MOS2 image in the three energy bands: soft (S; 0.2-1.0 keV), medium (M; 
1.0-2.0 keV), and hard (H; 2.0-12.0 keV), using the SAS task EDETECT_CH AIN . A likelihood 
value of 10 was used, corresponding to a significance level of 3.6a. The outputs from EDE- 
TECT_CHAIN were inspected for spurious sources at the edges and boundaries between chips. 
To generate the final source list, we imposed two additional selection criteria for sources detected 
by EDETECT_CHAIN for a source to be confirmed as a true detection if it satisfied either one of 
the followings: First, the source had to be detected in at least one additional observation within a 
searching radius of 6"; second, if the source was detected in only one observation, the S/N of the 
source had to be greater than 5 in that detection. These criteria were applied to all observations, 
except for the 2001 February ob servation in w hic h the sources found by EDETECT_CHAIN were 
consistent wi th the source list in lKongl(|2003|) and lBauer et al.1 (|2003|) and thus we just adopted the 
source list in lKongl(|2003|) for this observation. About one third of the sources detected in the four 
observations were eliminated by these criteria. 

Using this algorithm, we found 37, 43, 51, and 30 sources in the observation taken in 2001 
February, 2004 February, 2004 August, and 2005 February respectively. Combining the individual 
source lists, we identified a total of 61 sources in our final X-ray catalog. As a cross check, we 
compared the sources detected b y the EDETECT CH AIN software in SAS with those detected by 
the WAVDETECT task of CIAO dFreeman et alJ2002h . running on a merged MOS+PN image with 
a significance threshold of 1 x 10~ 6 . Most sources that were detected with EDETECT_CHAIN in 
the final catalog could also be detected by WAVDETECT, except for sources 1, 2, and 53 which 
were only detected with EDETECT_CHAIN in the longest exposure data (2004 Auguest). It is 
possible that they are spurious sources due to peculiar point spread function features or the presence 
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out-of-time event signatures (|Argiroffi et al.1 120061) . On the other hand, these sources might be 
transients that could be detected during limited period. However, the fact that the three sources are 
all observed only in the deepest observations could simply mean that they were below detection 
limit in the other exposures. Therefore we consider them as real detections. 

The absolute astrometry of source positions are limited by two factors. Firstly, the aspect so- 
lution of XMM-Newton is accurate to w 4". Second, the systematic shifts among the observations 
made it possible that source detection processes run in different bands and observations might give 
slightly offset centroids for the same source. In order to determine source positions, we first cre- 
ated a single- stacked image by aligning the four observations in the coordinate frame of the 2004 
February image. Eight bright X-ray sources common to all data sets within 10' of the aimpoint were 
selected for the alignment purpose. The relative shifts between these exposures were determined 
using the IRAF task CCMAR The cross registration gave a rms accuracy of 0.8". Positions of the 
61 discrete sources in our final X-ray catalog were then determined in the aligned single-stacked 
image with WAVDETECT. The positional errors from WAVDETECT varied as the off- axis angle 
from the center of the image, with sources at the outer edges having larger errors up to 2" and these 
were tabulated in Table [21 We attemped to improve the astrometic accuracy b y cross-correlating 
the output X-ray pos itions from WAV DETECT with the optical USNO-B1.0 JlVIonet et alJl2003h 
and near IR 2MASS (|Cutri et al.ll2003 ) catalogs. We identified 25 X-ray sources with optical coun- 
terparts and 16 having near IR counterparts, using a searching radius of 5". We selected 7 X-ray 
sources which possessed consistent and high S/N optical counterparts in at least three observa- 
tions, to correct for the astrometry of the X-ray images with the optical positions using the IRAF 
CCMAP task. The average astrometic rms were 1 .05" and 1.14" in RA and DEC respectively. The 
positions listed in Tableland throughout this paper are based on this astrometric reference. 

The spatial distribution of the X-ray sources in IC 342 is not uniform. The majority of the 
sources detected are located on the spiral arms of the gal axy (Figure [I])- The 17 sources that are 
outside theZ>25 disk (D ~ 21.4', de Vaucouleurs et all 1992 ) of IC342 are most likely to be contam- 
inated by foreground or background objects. To esti mate the number of b ackground objects within 
the D25 disk, we used the Chandra Deep Field data (|Brandt et al.ll200l|) and concluded that there 
were about 10 background sources based on a detection limit of fx ~ 2.7 x 10~ 14 erg cm -2 s -1 , 
corresponding to 23% of the total number o f sources inside the D25 disk. This number is con- 
sistent with the result by iBauer et al.1 (|2003|) (7-12 ba sed on the detection limit of fx ~ 1.0 x 



10 14 erg cm 2 s l ) but higher than that by llCong ( 2003 ) (~ 3 based on the completeness limit of 



fx ~ 7.7 x 10 14 erg cm 2 s 1 ). The difference in the number of background sources estimated 
from these two groups is probably due to the use of different flux limit values. 
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3.2. Cross-correlation with X-ray catalogs 



We performed a systematic cross-correlation with existing IC342 catalogs from other X-ray 
missions including Einstein, ROSAT, Chandra, as well as previous XMM-Newton publications. 
Source identification are compared with this new catalog by varying the searching radius according 
to the spatial resolution of the corresponding X-ray instruments. The results were listed in Table [3] 
The 3 brightest X-ray sources in IC342 ULX-1, ULX-2, ULX-3, first detected by Einstein, were 
also detected by t he current XMM-New ton observations. For the ROSAT sources, all except R2 
(designation from lBregman et al.lll993|) were detected in our current XMM-Newton observations. 
R2 was the faintest source in the ROSAT observation with only 14.7 ±5.5 counts in the 0.1-2.5 keV 
range, suggesting that it might be too faint to be detected in our observations and was thus excluded 
in our source detection algorith m. We found 21 of our XMM-Newton sources to have counterparts 
in the Chandra HRC-I catalog (IMak et al. hood) . Four of the Chandra sources were not detected in 
our XMM-Newton observations: two (CIO, C13; designation from LMak et al.ll2008|) of them were 
located close to the nucleus and were not resolved in the XMM-Newton data; one (C20, probably 
a foreground star) of them was outside the field of view; and the remaining one (CI 9) was either 
below the detection limit of XMM-N ewton or might be a varia ble. We also recovered all the XMM- 
Newton sources previously listed in iBauer et al. 120031) and Kongl (|2003[) that was based on the 
2001 February observation. The sources lis t s in b oth works were very similar, with iKongl (|2003|) 
listing two more sources than IBauer et al. I J2003L The X-ray properties, including color-color 
diagram, spectral fits of bright sources, of the 2001 February observation analyzed in our work 
were consistent with these two references. 



There were 22 XMM-Newton sources that were not detected in previous Einstein, ROSAT, 
and even Chandra observations, presumably because they were below the detection limit or not 
resolved by these instruments. It is worth noting that the XMM-Newton observations in the present 
analysis should reach the lowest flux level. For example, the exposure time of the 2004 August 
XMM-Newton observation at 23.6 ks was a factor of ~ 2 times longer than the Chandra HRC-I 
observation, and 20% longer than the ROSAT observation. Nevertheless, it is possible that these 22 
sources might be variables and were thus not observed in previous X-ray data. As we will discuss 
later in section |4~T1 14 of these 22 sources were classified as variables and/or X-ray transients. 
For example, source 45 was detected in both the 2004 observations and became very luminous, 
with Lx = 4.9 x 10 38 ergs _1 , but it was dimmer than the detection limit in the 2005 February 
observation. These comparisons helped to identify the long-term variability of the X-ray sources. 
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3.3. Photometry 



There are two methods to perform photometry on the X-ray images: one with the SASDAS 
task EMLDETECT specifically designed for analyzing the XMM-Newton data and another one 
with the CIAO task DMEXTRACT that is compatible with both XMM-Newton and Chandra data. 
The count rates and fluxes from the output of EMLDETECT have been corrected for vignetting 
and instrumental PSF that account for the variation across the CCD before they are background 
subtracted. However, there is a limitation in using this tool to determine the net source counts 
since we do not have detailed knowledge of the background subtraction and exposure correction 
implemented within the tool (M. Ehle, private communication). Therefore some of these source 
count rates could possibly be contaminated by nearby sources in crowded regions (e.g. the center 
of IC342). On the other hand, the net source count generated from the task DMEXTRACT do not 
suffer from this problem. As a result, the net source counts were determined by DMEXTRACT 
while the count rates and fluxes were determined by EMLDETECT in this work. 

The X-ray fluxes of individual sources were determined in each observation, even when the 
source was below the detection limit in that particular data set. This was done by running the 
SASDAS task EMLDETECT with the final position of each source as an input parameter (by set- 
ting fitposition=no) and setting the detection threshold zero. The aperture sizes were automatically 
determined by the maximum likelihood PSF fitting algorithm. Furthermore, to convert the count 
rates to fluxes, we used the energy conversion factors (ECF) calculated using WEBPIMMS by 
assuming an absorbed power-law spectrum with a photon index of 2 and an absorption column 



density of 8 x 10 21 cm~ 2 (typical for X-ray Binaries: lKongll2003|) . The resulting ECF values were 
listed in Table H] and were used as input parameters for EMLDETECT. This power-law model was 
used since the X-ray population is dominated by X-ray Binaries and the fact that hardness ratios 
of most sources in IC342 are consistent with such model (c.f. Section |3~4l ). For non-detections, we 
determined the 3a upper limit of the flux and luminosity of the sources assuming a distance of 1.8 
Mpc to IC342. 

To extract source counts, we produced mosaiced X-ray images from each observations that 
combined data from the three EPIC cameras in the full band, followed by normalization using the 
exposure maps created from the SAS tool EEXPMAP. We extracted the source counts via aperture 
photometry and varied the aperture size for each sources based on their off-axis angle in order 
to match the 90% encircled energy function. The extraction radii were smaller near the aimpoint 
where the PSFs were well defined and bigger at larger off-axis angle where the spatial resolution 
was poorer. The net counts for each of the three bands listed in Table |2] are the sum from the 4 
observations. 
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3.4. Hardness Ratio 



Most of the X-ray sources in IC342 are faint with fewer than 200 counts in each observation, 
which makes it difficult to derive accurate spectral parameters. Hardness ratios can provide crude 
indications of the X-ray spectra in these cases. We computed the hardness ratios for each detected 
source, based on the sourc e counts in soft (5), medium (M), and hard (H) bands. Following the 



convention in 



Kongl (|2003l) . we computed two hardness ratios defined as HR\ = (M-S) / (M+S) 
and HR2 = (H-S)/(H + S). Table [2] lists the average values for both HR1 and HR2 over the 
four observations for each detected sources. Figure [2] shows the color-color diagram (CD) for 
all detected sources. We also checked the CD for each of the four observations individually and 
fo und that they were similar. In particular, the CD from the 2001 data is consistent with the results 



of iKong d2003k We overlaid on the CD six lines showing the tracks followed by representative 



spectral models with different values of N^: power-law models with power-law index T = 1.2, 1.7, 
2, and 3; Raymond-Smith model with kT^s = 0.5 keV; blackbody model with kT^B =0.1 keV. 
For each model, A^h varied from left to right with 3 x 10 21 cm~ 2 , 5 x 10 21 cm~ 2 , and 10 22 cm~ 2 . 
The power-law spectra tend to occupy the upper right region of the CD while the soft thermal 
models tend to occupy the lower left region. Typical X-ray binaries or AGNs with power-law 
spectrum would most likely be located at the upper right region while supersoft sources (SSS) 
having blackbody spectra of kTsB ~ 0. 1 keV would occupy at the lower left. 

The X-ray sources in IC342 are dominated by sources emitting power-law spectra, with only 
a few of them compatible with thermal spectra. These sources with soft X-ray emission could pos- 
sibly be supersoft X-ray sources, supernova remnants, or foreground stars. This trend is consistent 
in all the four datasets as shown in their individual CDs. To further classify the XMM-Newton 
sources into subcategories s uch as foreground stars , SNR, and SSS, we employed the hardness 



ratio selection a lgorithm by iMisanoyic et al.l ((2006) for the foreground star and SNR classifica 



tion, and that of Di Stefano & Kongl l2003bl) for SSS and Quasisoft sources (QSS) classification. 



Since the definition of the hardness ratios employed in these publications are slightly different from 
our chosen energy bands, we computed the corresponding HR values used for each classification 
scheme. Detailed results will be discussed below in Section |4~3l 



4. Results 



4.1. Flux Variability 



The four XMM-Newton observations analysed in this work spanned a period of four years from 
2001 to 2005, allowing us to study the long-term X-ray variability of these sources. Following the 
definitions by IPrimini. Forman. & Jonesl (|1993|) . we computed a flux variability parameter 5fl ux 
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defined as 



l^rnax ^min 



ri | iiioa 1111111 / 1 \ 

^flux = 1 = , U) 

\ (7% + o\ 
y r mas ^min 

where F max and F m j n are the maximum and minimum X-ray full-band fluxes of the four obser- 
vations with <jp , and ok as their corresponding errors. A source is defined to be a flux 

r max r min A ° 

variable if 5fl ux > 3 and is marked with "f" in the last column of Table [2l If the minimum flux 
is a non-detection (i.e. S/N< 5), we used the 3a flux upper limit to calculate the lower limit 
of the variability factor. Since most of the sources in our catalog have low number of counts 
with nonzero background, we empl oyed the Bayesian approach to calculate the 3a flux upper 



limit (Kr aft. Burrows. & Nouseklll99ll) . 



Long-term X-ray flux variability were detected for 39 out of 61 sources , representing 64% 
of the total. Of these, 29 sources varied in flux by more than a factor of two (i.e. F^x/F^ > 2). 
The maximum amount of variabiliy corresponds to 5fl ux ~ 60. We plotted the relation between the 
average offset from the galactic center and the variability factor of each source in Figure [3] It is 
noted that most of the strong variables are located between ~ 3' -1' , corresponding roughly to the 
region of the spiral arms, hinting that these sources are indeed located within the galaxy and are 
not foreground or background sources. 



We have also adopted the criteria of iKong et al.1 (I2002b|) to search for X-ray transients which 



are defined as sources: (1) having flux variability factor 5fl ux > 3, and (2) found in at least one 
observation with Lx > 10 37 ergs _1 and was not detected (i.e. source counts are below the 3a 
detection threshold) in at least one other observation. The luminosity limit covers typical outburst 
luminosities of soft X-ray transients and Be/X-ray binaries in our Galaxy. A total of 26 X-ray 
transients are detected (about 43% of the total 61) and are marked with "t" in Tabled About half 
of these transients were bright during their high state, Lx > 3 x 10 37 erg s -1 , the brightest of which 
(source 36) has Lx = 4.9 x 10 38 ergs _1 , meaning that the source had brightened by a factor of 
> 50. 

We plotted the long-term 0.5- 10 keV light curves for the three ULXs (sources 19, 25, 38) 
and one bright source (source 17) in Figure HI A distance to IC342 of 1.8 Mpc was assumed and 
the luminosities were not corrected for absorption along the line of sight. The baselines of the 
light curves were extended by adding data points from the 1998 ROSAT HRI results and the 2006 
Chandra HRC-I observations. These data were extracted from the archive and were reanalyzed 
using standard pipeline procedures of the corresponding instruments. The derived luminosities in 
this work cannot be compared directly with the previously published ones because of the different 
assumptions made in those works. For example, luminosity values in the earlier ASCA and ROSAT 
publications were computed assuming a distance to IC 342 of 4.5 Mpc instead of the 1.8 Mpc 
assumed here. We recalculated luminosities of all sources from all observations by assuming 
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the same energy range, spectral model, and assumed distance of IC 342. Identical extraction 
aperture for each source had been used in the XMM-Newton, ROSAT, and Chandra data, which 
was reliable for source 17, 19, and 25. However, with source 38 located in the crowded region 
near the center of IC 342, it suffers from contamination of nearby sources as well as the diffuse 
background of the center. Using a r = 15" circular aperture centered on source 38, it was estimated 
that roughly 40% of source counts were contaminated by the faint source C13 (IMak et al.l 120081) 
which was unresolved in o ur XMM-Newton i mages, and the local background when compared to 
the photometry obtained by lMak et al. This suggests that the flux values estimated for this 

source from XMM-Newton and ROSAT data are overestimated due to their larger PSF and as a 
result we should keep that in mind while examining its light curve. 

The flux values plotted in Figure 0] were obtained from spectral fits to the four observations 
using XSPEC. For consistency, the fluxes of the four XMM-Newton observations of each sources 
were derived from the same spectral model with identical best-fitting parameters (Table [8]). On 
the other hand, with no spectral information in the ROSAT HRI and Chandra HRC-I observations, 
the fluxes were instead estimated using PIMMS using spectral parameters averaged from the four 
XMM-Newton observations. With the exception of source 38 which is known to be a spectral vari- 
able (see Section 1431) . the spectral parameters of sources 17, 19, and 25 do not vary much among 
the four observations and thus the uncertainties in the derived Chandra and ROSAT fluxes from 
these parameters for each source is small or negligible. For source 38, the Chandra and ROSAT 
fluxes derived from the four sets of spectral parameters (using the power-law + blackbody + gaus- 
sian model) vary by a factor of ~ 1.3. The implied flux ranges were accounted for in the error bars 
while the average value was plotted as the data point of its long-term light curve. For the other 
three sources in Figure H the error bars in the light curves for Chandra and ROSAT data points 
are dominated by poisson errors in the counts, while those in XMM-Newton data points are la 
gaussian confidence limits from XSPEC. It should also be noted that comparisons of luminosities 
between different instruments should be treated with caution because various factors (e.g. cross- 
calibration issues and uncertainties in instrument responses and in assumed spectral shapes for 
energy conversions) could possibly contribute to the errors. Previous studies indicated t hat varia 



tions of normalized flux in different X-ray instruments were small, at a level of ±10% (|Snowden 



20021), an d thus were not expected to account for all the flux variability observed. 



4.2. X-ray Spectral Analysis and Spectral Variability 



We searched for spectral variability of our XMM-Newton sources using two methods. We first 
inspected changes in the best-fit spectral shapes and/or measured changes in the best- fit parameters. 
Similar measurements of spectral variability had been used by iGrimm et al.l (|2007[) to search for 
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spectra variables in M33. However, this could only be done for bright sources that possessed 
enough source counts for reliable spectral fits. Second, we adopted a spectral variability factor 
analogous to the flux variability factor described in Section [47TI with flux in equation (0Q) replaced 
by the hardness ratios. A source is defined to be a spectral variable if 5 spec trai > 3. This method 
could be performed on all sources, including the faintest ones. 



4.2.1. Spectral Fitting 

We extracted the energy spectra of the 1 1 brightest sources from all our four XMM-Newton 
catalog with the SAS task XMMSELECT. These 11 sources included those that were highly 
variable in flux or spectral. Response matrices were created by RMFGEN and ARFGEN while 
X 2 statistics was used to find the best-fitting parameters. We first tried simple absorbed one- 
compon ent models including power-law, blackbody, multicolor disk black-body (diskbb model in 
XSPEC; lMitsudaetal.1l 19841) . Raymond-Smith (RS), BMcfl, and broken power-law models. For 



a few cases two-component models were attempted when single-component model did not gener- 
ate satisfactory fits. Spectra from the three XMM-Newton detectors were fitted simultaneously for 
each source, allowing only an overall normalization factor to vary among detectors to account for 
differences in calibration. 

The spectra of these 11 sources were satisfactorily fitted with simple absorbed power-law 
models, except for source 27 (SNR candidate), 36, and 38 (the nuclear source). The best fit spectral 
parameters are listed in Tabled The best fit Nn ranges from 2.3 x 10 21 cm~ 2 to 2.5 x 10 22 cm~ 2 , 
with an average value of 8.0 x 10 21 cm~ 2 , while the photon index varies between 0.47 and 3.52. In 
general, at least 1000 detected counts were needed to clearly rule out competing spectral models. 
In spite of this, the deduced fluxes do not vary significantly among different spectral models and 
thus it is possible to derive X-ray luminosities and study long-term flux variabilities. Detailed 
fitting result of three brightest X-ray sources (sources 19, 25, 38) and the SNR candidate (source 
27) will be discussed below in detail in section I5T21 

Five sources (19, 23, 24, 25, 38) were classified to be spectral variables based on spectral 
fits. For sources 19 and 25, their spectra shapes changed significantly not only among the four 
XMM-Newton observations, but also over the past 10 years when compared to previous ASCA 
observations. The changes are reminiscent of spectral transitions of black hole binaries (BHBs) 
which usually vary between the high/soft and the low/hard state. For source 38, the spectrum in 



3 The typical scenario involved thermal X-rays from the inner region of an accretion disk in a 
black-hole binary illuminating in-falling matter in close proximity to the black-hole event horizon. 
(http://heasarc.nasa.gov/xanadu/xspec/manual/XSmodelBmc.html) 
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2001 was significantly different from the other three with a larger photon index at energy above 
1 keV. We shall discuss in detail the X-ray properties of these three sources in Section [5721 The 
fitted spectral parameters of these five sources all showed variability at the 90% confidence level 
(corresponding to a 2.106a confidence limit). For example, the power-law fits of source 24 for the 
2004 February and 2004 August data resulted in a change of photon index Ar = 2.4, implying a 
photon index change at more than 99% confidence level (corresponding to a 3a confidence limit). 
The power-law + MEKAL fits to data of source 38 also showed a trend of decreasing photon index 
(Ar = 1.3) and increasing plasma temperature (AkT « 0.37) at the 90% confidence level. 



4.2.2. Spectral Variability Factor S S pectrai 

We define two spectral variability factors Sspectral by replacing the flux with HR1 and HR2 in 
equation ([]])• High spectral variability with S S pectrai > 3 in either HR1 and HR2 was detected only 
in three of the sources (source 19, 36, 41) in which two (19 and 36) had already been classified 
as spectral variables from spectral fits. An alternative method to search for spectral variable is to 
study the color- luminosity diagrams. We inspected the color-luminosity diagrams of all sources 
and found that 15 of them exhibited significant spectral transitions in different patterns. Five of 
them (source 20, 32, 40, 45, and 51) showed typical spectral behavior reminiscent of Galactic X- 
ray BHBs which were seen as having negative slopes in the plots; one of them (source 52) showed 
correlation between lower luminosity and decreasing hardness; eight of them (source 9, 16, 27, 
30, 34, 40, 43, 55, 56) showed more complicated behavior or underwent several transitions; one 
of them (source 38) varied spectrally but not in intensity. As a result, a total of 19 sources, or 
~ 30%, were classified as spectral variables under these criteria and they were mar ked with "s" 



in Ta ble [2] Spectral variability seen in X-ray sources is common in external galaxies. IZezas et al. 



(120061) found 21 out of 120 X-ray sources, ~ 18%, in the Antennae show signs of variability in 



hardness ratios while IFridriksson et al.1 (|2008) found that 2 out of 90 and 3 out of 38 X-ray sources 
in NGC 6946 and NGC 4485/4490 respectively show significant (3cr)variability in the hardness 
ratios as defined similarly in this work. It is also noted that while these galaxies, just like IC342, 
are starburst galaxies, the fraction of spectral variables we identified in IC342 is higher than these 
galaxies. 



4.3. Identifications of X-ray Sources 

Using properties of the detected sources (e.g. HR, variability) and cross-correlating them with 
catalogs at other wavelengths, we attempted to classify all sources in our X-ray catalog into dif- 
ferent classes of X-ray emitting objects. We varied the searching radius according to the accuracy 
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of the various catalogs and visual inspection was performed to confirm the spatial coincidence of 
a counterpart. A summary of our classification scheme is presented in Tableland the results of 
the source identification are listed in Table [7J We briefly explain the classification of each type of 
objects below. 



4.3.1. Super soft sources and Quasisoft sources 



Supersoft X-ray sources (SSS) are characterized by their soft X-ray emissions, which can be 
represented by blackbody spectrum with kT < 100 eV, and with X-ray luminosities of 10 36-38 erg s -1 



(|Di Stefano & Kongll2003a|) . These sources show little or no X-ray emission above 1-2 keV, and 
are generally believe d to be binary systems with white dwarfs accreting from more massive hydro- 
gen burning donors (Ivan den Heuvel et al.N 19920 . Some of them show significant v ariabilities on 



vario u s time scales, while some other have been found to be recurrent trans ients (|Greiner et al. 



20041 : iDi Stefano et al.l 120041 : iKong & Pi Stefanol 120031 : lOsborne et al.l 1200 II) . Quasisoft X-ray 

sources (QSSs) are luminous (L x > 10 36 ergs _1 , kT between 120 eV and 350 eV) X -ray sources 

emitting few or no photons at ener gy above 2 keV yet clearly emitting at above 1 . 1 keV (IDi Stefano & Kong 
2003bl : Pi Stefano & Kong Il2004|) . They also suggested that if we observe a hot SSS located be- 
hind a large gas column, just like the case of Galactic center, photons in the medium energy band 
(1.1-2 keV) would be detected, while only few soft photons would be detected. These sources 
might be the hottest nuclear burning white dwarf binaries and could possibly be progenitors of 
Type la supernovae, and SNRs. 

Since most sources in our catalog had too few photons (counts < 200) f or meaningful spec- 
tral fitti ng, we employed the selection algorithm for SSS and QSS defined by IDi Stefano & Kong 
(I2003bl fig. 3) based mainly on hardness ratios. The energy bands defined in the hardness ratios of 
their work were different from what we used earlier in Table |5l in which the three energy bands 
were defined as: soft (5) =0.2-1.1 keV, medium (M) = 1.1-2.0 keV, and hard (H) = 2.0-7.0 
keV. We therefore recalculated hardness ratios according to these definitions. The summary of 
source counts and hardness ratios of the identified SSS and QSS is given in Tabled 

Excluding those that had been classified as foreground stars or candidates, we identified 
source 36 and source 56 as confirmed SSS and QSS respectively, and we describe source 36 im- 
mediately below in detail. Source 56 satisfied the QSS selection criteria in the 2002 February data 
but was too faint to be identified as a source in the remaining three observations. Two additional 
sources (52, 54) also satisfied the QSS selection criteria. Source 52 satisfied the selection criteria 
in 2005 February in which the S/N is the lowest among the four observations. It however pos- 
sessed substantial hard signal in the other three observations. Source 54 was identified based on 
the 2004 February and 2005 February observations in which the X-ray emission are dominated by 



13 



photons below 2 keV. On the other hand, substantial hard signal was present when the source was 
in its high state in 2004 August. We therefore classified both sources 52 and 54 as possible QSS 
candidates and also presented their properties in Tabled Interestingly, all the classified QSS and 
QSS candidates are under the category of "QSS-noh", meaning little or no hard X-ray emission 
was detected. 

The only classified SSS, source 36, had also been identified as a transient. It was very lu- 
minous during its high state in both observations in 2004 but disappeared from the field in the 
2005 February observation, just after it reached the highest state in 2004 August. Spectra of 
the source were extracted from the two 2004 observations. Simple models all gave unaccept- 
able fits with x 2 > 2. Satisfactory fits could be obtained only after a gaussian line at « 0.9 keV 
was added to the blackbody and power-law model. A relatively high temperature of kT\^ 0.17 
keV was obtained for the BB model, while very large photon indices of T — 3.6 and 6.8 were 
obtained for the power-law model with the 2004 February and August data respectively, indicat- 
ing soft X-ray emissions. We noticed that SSS transie nts are common in nearby galaxies, with 
two luminous examples being source 1 10 in N GC4697 dSivakoff et al.ll2008|) and ULX-1 in M101 



(IKong etal.ll2004al : 



Kong & Di Stefano Il2005|) . They both belong to the group of ultraluminous 



SSS (ULS) that have been suggested to be accreting intermediate mass black holes (IMBHs) 
(IDi Stefano & KongH2003a|) . In our case, source 36 is likely to be consistent with the white dwarf 
model since its luminosity conforms to the near-Eddington luminosity of a 1.4M Q white dwarf 
(L^ dd ~ 1.8 x 10 38 erg s -1 ). We looked for its optical counterpart from the Digitized Sky Survey 
(DSS) data and UV/optical counterpart from the XMM-Newton Optical Monitor (OM) images, but 
found no source at its position in all images. 



4.3.2. Supernova Remnants 



With IC342 not being an active target for SNR survey, there had been few identified SNR 
candidates in the galaxy. We searched for SNR and stellar novae as listed in the International 
Astronomical Union Circulars (IAUCs), and no matches within 10' of IC 342 was found. On 



the o ther hand, we found four SNR candidates listed in the optical search by iDodorico et al 



(I1980|) . The SNR object 1 (based on designation of IDodorico et al.l 1 1980) is 9.4" away from 
our X-ray source 27, with uncerta inty in position of th e SNR at about 10" from the optical 
data, while the other three SNR in IDodorico et al. Jl98Ch were not within the proximity of any 
of our XMM-Newton sources. This SNR was first identified in radio continuum observations of 



the galaxy (Ba ker et al.lll977l). and was later confirmed as a SNR on the basis of optical spec- 
troscopy ( Dodorico et al.l 1980 ). Detailed X-ray analysis of this source will be presented in Sec- 
tion !5.2.4l Its X-ray spectrum could be fitted with an absorbed NEI model (c.f. Section 1531) with a 
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derived 0.5-10 keV luminosity at ~ 9 x 10 36 erg s" 1 , which is similar to that of SNRs in the Mag- 
ellanic Clouds (e.g. iHughes. Hayashi. & Koyamall 19981 . Iwilliams et al.ll 19990 . Its X-ray spectrum, 
combined with its proximity to an optical SNR suggest that source 27 is a strong SNR candidate. 

In addition to comparison with o ptical SNR catalog, we also attempted to identify SNR using 
the hardness ratio criteria outlined in lMisanovic et al.1 (120061) . which required a SNR candidate to 



have//#l z > 0.1 and//i?2 z < -0.4, where HRl z = (S 2 - Si) / {S 2 + Si) mdHR2 z = (5 3 -5 2 )/(5 3 + 
^2), with Si, S2, and S3 the source counts in the energy bands 0.2-0.5 keV, 0.5- 1.0 keV, and 
1.0-2.0 keV respectively. Excluding the foreground candidates (see below), two sources (48, 59) 
were found to satisfy this criteria. It was noted that the aforementioned SNR candidate source 27 
did not satisfy the HR criteria. As we will discuss in Section I5.2.4L the X-ray spectrum of source 
27 could also be satisfactorily fitted by a power-law model and thus it was possible that a hard 
component resided in i ts X-ray emission. Together with the fact that the hardness ratio criteria of 
Misanovic et al.1 (|2006|) only considered X-ray emission below 2 keV, this source could therefore 
possibly be missed. It also implied that similar sources with a relatively hard component could be 
missed with this criteria as well. 

For the two SNR candidates selected by the HR criteria, we checked for their variabilities 
as outlined in Section |4T1 We found that source 59 exhibited some degrees of variability in flux, 
with Sfl ux = 4.7. Since the X-ray emission of a SNR is expected to be persistent, this source was 
rejected from being a SNR candidate. On the other hand, source 48 had a very soft X-ray emission 
which was consistent with a QSS classification. It had an optical counterpart in the USNOB1.0 
catalog with an offset of 0.58". Moreover, the X-ray to optical flux ratio at log(/ x // op t) = -0.3 
was larger than the exp ected value for foreground stars, which should have log(/ x // pt) < _ 1 
(IMaccacaro et al.lll988|) . Since optical emission of accreting binary systems at the distance of 



IC342 is below the sensitivity of the USNO and 2MASS catalogs, therefore source 48 was classi- 
fied to be a SNR. 



4. 3. 3. Foreground Stars 

We compared our XMM-Newton source list with the USNO B1.0 catalog and found that 25 
of them had optical counterparts within a 5" searching radius (except for source 59 which had an 
offset of 6.2" with an USNO star, but the searching radius was relaxed because it was at a large 
off-axis position). We found 13 of these 25 sources have X-ray to optical flux ratios that were 



consis tent with those of normal stars of log(/ x // op t) < 1, using the criteria of IMaccacaro et al. 



(119880 . We calculated the X-ray to optical flux ratios from log(/ x // opt ) = log(/x) +0.4V + 5.37 for 
each source. In the calculations, the X-ray flux f x was calculated by assuming a simple power-law 
model restricted in the energy range 0.3-3.5 keV with T = 2 and Nu = 8 x 10 21 cm~ 2 (same model 
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as assumed in luminosity values in Table [2]), and the value of the V magnitude was averaged from 
the B and R magnitude in the USNO B1.0 catalog. To further constrain the classification of these 
13 sources as foreground stars, two additional criteria are proposed: 



1. Hardness Ratio: X-ray emission from stars are relatively soft from studies of their spec- 
tra or hardness ratios. The energy spectrum of a foreground star could be best fitted by 
a Ray mond- Smith model with &7rs ~ 1 keV or a power-law model with a photon index 



r > 3 (kong et al.ll2002b|) . In our XMM-Newton catalog, 13 sources were n ot bright enough 



for sp ectral analysis. Instead we used the hardness ratio criteria defined in lMisanovic et al. 



(2006), which required a foreground stars to have HR2 Z < 0.3 and HR3 Z < -0.4, where 
HR3 Z = (S4 — S3)/ (S4 + S3) and 54 was the source counts in the energy band 2.0-4.5 keV, 
while S3 and HR2 Z were the same as that defined in Section l4.3.2l 

2. 2MASS counterpart: For sources that had near-IR counterparts from the 2MAS S catalog, we 



check ed for their near-infrared colors and magnitude. Adopting the criteria of iFinlator et al 



(I2000|) . we classified sources with J—K < 0.8 and J < 12.5 as foreground stars. 



To classify sources as foreground stars, four criteria had been proposed in which the first 
two, that is, USNO counterparts and X-ray to optical flux ratios, were used for screening purpose. 
Based on the two additional criteria, we classified X-ray source to be a foreground star of category 
1 ("Catl " in Table[7]) if all four criteria were fulfilled. Sources that satisfied either the hardness ratio 
or the 2MASS criteria in addition to the two screening criteria were identified as foreground star 
candidates and were marked as categories 2 and 3 respectively ("Cat2" and "Cat3" in Table©. We 
confirmed four sources (7, 24, 26, 29) as category 1 foreground stars that sati sfied all the criteria 



described above. These four sources were also identified as foreground stars by lBauer et al.1 ((2003) 
based solely on the X-ray to optical flux ratio. Furthermore, we identified seven other sources for 
the first time as foreground star candidates since they satisfy either the hardness or the 2MASS 
criteria. The majority of these candidates, six in total (source 16, 28, 37, 43, 51, 61), belonged to 
the category 2, with only one (source 2) belonging to category 3. 

In summary, a total of 11 sources were classified as foreground stars or candidates. They all 
had bright optical counterparts, with visual magnitude B < 20.2 and R < 17.1. Except for sources 
2, 16, and 61, they all showed significant X-ray variability on the time scale of years, probably due 
to flarings. In particular, sources 43 and 51 were transients and exhibited large changes in flux. 
Source 24 was the only star candidate that had enough counts for reliable spectral fitting. Both 
the Raymond-Smith and power-law models could generate satisfactory fits to the spectrum, with 
&7rs ~ 0.9 and power-law index > 3 (c.f. Table [8]), which was consistent with the spectrum of 
foreground stars. 
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5. PROPERTIES OF X-RAY SOURCES 



5.1. Global Properties 



Analysis of the four observations of IC 342 with XMM-Newton between 2001 February and 
2005 February were presented in previous sections. The limiting luminosity of these exposures was 
between ~ (1 -2) x 10 37 erg s -1 depending on the local background and exposure time. Combining 
all four observations, a total of 61 X-ray sources were detected. While most detected sources 
were too faint for detailed analysis, we used the X-ray colors, flux and spectral variability, and in 
some cases optical or radio counterparts to classify the detected X-ray sources into several stellar 
populations. We identified a total of 1 1 foreground stars, 1 SSS, 3 QSS (one confirmed and two 
candidates), and 2 SNRs, with a summary given in Table[7J The spatial distribution of all detected 
sources was shown in Figured! with classified sources characterized in different symbols. The 
SSS/QSS were found to be located near or on the spiral arms, associating them with young stellar 
populations. On the other hand, the two SNR are very close to the starburst nucleus dominated by 
current star formation activities. 

One major focus of this work is the study of intensity and spectral variability of the X- 
ray sources in IC342 on timescales of years. We found that 41 of the 61 detected sources, or 
64%, showed long-term flux variability, clearly indicating that they were individual X-ray bina- 
ries. Of these, 26 sources, or ps 43%, were classified as X-ray transients. The observed fraction 
of sources showing flux variability of (64%) is quite high compared to previous studies of late- 
type galaxie s. For example, 27 % of the X-ray sources in the 11 nearby face-on spiral galaxies 
studied by iKilgard et al.l (|2005r) exhibited variability on either long or s hort timescales, 25% of 
the sources in M33 exhibited long-term variability (IGrimm et al.1 12007|) . 29% of the sources in 
NGC694 6, and 39% of the sour ces in NGC4485/4490 were variables on timescales from weeks 
to years (|Fridriksson et al.1 120081) . This indicated that th e source pop ulation in IC342 was domi- 
nated by accreting XRBs, in agreement with the results of lKongl (|2003|) that sources in IC342 were 
mostly HMXBs. In addition to flux variability, 19 sources, or 30% of total, were identified to be 
varying in hardness ratios or undergoing spectral transitions. Eight of the identified X-ray tran- 
sients showed spectral variability in addition to flux variability as seen from the color-luminosity 
diagrams. 

In Figure [5l the X-ray hardness ratios (HR1 and HR2) versus the 0.5-10 keV X-ray lumi- 
nosities of all detected sources were plotted. Most sources were in the luminosity range 10 37 - 
10 38 erg s l , and their hardness spanned a diverse population. The HR1 color-luminosity diagram 
(upper panel) showed that a majority of the sources are hard while the HR2 diagram (lower panel) 
revealed a large proportion of hard sources along with a separate and smaller soft sources. This 
soft population in HR2 color comprised of foreground stars, SNRs, SSS/QSS, with only one source 
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that was unclassified. This is consist ent with the expec tation that these sources are characterized 
by their soft color. As pointed out in lZezas et al.1 (|2006r) . the observed HR1 color was sensitive to 
absorption while HR2 was most sensitive to the intrinsic spectral shape. With the X-ray sources 
in IC342 being mostly hard and slightly obscured, this picture is consistent with a population 
dominated by HMXBs, as would be expected in si tes of recent star formation. For a compari- 
son with an extragalactic X-ray source population, llrwin. Athey. & Bregmanl (12003b studied 15 
nearby early-type galaxies observed with Chandra and found that sources with luminosities in the 
(1-2) x 10 39 ergs _1 range had softer spectra (power-law T ~ 2), which was consistent with the 
high/soft state of black hole binaries. On the other hand, the two sources in IC342 with luminosities 
above 10 39 erg s -1 were instead generally harder, with hardness ratios > 0.5. 



5.2. Individual Sources 



5.2.1. Source 19 (IC342 X-l ): An ultraluminous compact X-ray source 



Sour ce 19 (IC342 X-l) is the mo st studied X-ray source in IC342. It was discovered by 
Einstein (Fabbi ano & Trinchieril Il987|) and confirmed as a point source with ROSAT observa- 
tions (jBregman et al.lll993r) . Subsequent ASCA observation in 1993 September sho wed that this 
sourc e was in a very luminous state with absorbed Lq.5-io keV = 3.8 x 10 39 ergs _1 (jKubota et al 



20021) (scaled to distance 1.8 Mpc) and its X-ray spectrum was best-fitted by an absorbed multi- 
color disk model (disk blackbody model in XSPEC), with 7V H = (4.7 ± 3) x 10 21 cm~ 2 an d 
7i n = 1.77 ±0.05 keV, which is generally used to describe ULX spectra ( Okada et al.l Il998|) . 
This corresponded to the characteristics of a black hole accretion disk in the high/soft state. The 
source dimmed by a factor of three in a follow-up ASCA observation in 2000 February with ab- 
sorbed Lo.5-10 keV = 1.2 x 10 39 ergs _1 (scaled to distance 1.8 Mpc) and the spectrum changed 
dram atically to an absorb ed power-law model, with Nn — (6.4 ±0.7) x 10 21 cm~ 2 and T = 1.73 ± 
0.06 (|Kubota et al.l l2002f) . The power-law-spectra as seen in many ULXs can be explained by 



the l ow/hard state that are observed in many Galactic and Magellanic BHB systems. LKubota et al 



(120021) noticed that there was significant softening in the X-ray spectrum above 5 keV and instead 
proposed that it was associated with an anomalous very high (V H) state (or recently described as 



steep power-law (SPL) state b ylRemillard & McClintockll2006l) . als o seen in many Gala ctic black 

l.hvlar 



hole binaries (e.g. GX 339-4. Markoff et al.ll2003l : GRS 1915+105. belloni et al.lboooh . The VH 
state is characterized by strong Comptonization and the ASCA spectra in 2000 could indeed be ad- 
equately fitted by a strongly Comptonized optically thick accretion disk with 7] n = 1.1 ±0.3 keV 
and the exponent of the radial dependence of the disk temperature r tn = 2.2 ±0.4. 



Using the XMM-Newton data from 2001 February to 2005 February, the 0.3-10 keV spectra 
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of source 19 were fit with absorbed power-law models (Figure [61 top-left) with spectra parameters 
listed in Table [8] The model gave satisfactory fits to all these observations except for the longest 
exposure 2004 August data. The photon index and the absorption column density increased from 
2001 February to 2004 February, followed by a slight drop six month later in 2004 August, while 
the inner disk temperature showed the opposite transition during the same period. It is worth to 
note that the spectral parameters of the disk blackbody model fit in the 2005 February data are 
remarkably similar to that of the 1993 ASCA observation. The spectral parameters from the power- 
law model are also consistent with those obtained from the 2000 ASCA observation. However, if we 
restricted the analysis to the energy range 2-10 keV, the power-law fits to the data showed that only 
the photon index in the 2001 February data was steep enough (r « 2.4), to be marginally consistent 
with the spectrum of t he VH state as defined by the pre sence of power-law component with T > 2.4 
in the 2-20 keV band ( Remillard & McClintocklEoOfjl) . The photon indices of the subsequent three 
observations were all lower than 2 in the 2-10 keV band. For the X-ray luminosity, the values 
in the energy range 0.5-10 keV also changed dramatically over the ~15 years. The long-term 
lightcurve of IC342 source 19 in Figure |4] (top-right) showed that it was one of the most variable 
sources in IC342 and its luminosity i ncreased by a f actor of more than 3 from 1991 to 2006. 
We noted that (and also pointed out by lKongll2003l and Bauer et al.ll2003l) the ASCA observations 
suff ered from serious c onfusion problem. The large extraction radius (3') used in analyzing source 
19 (|Okada et al.lll998|) would have included nine other sources (source 17, 18, 21, 24, 26, 27, 
28, 29, and 30) detected in our four XMM-Newton observations. Using the present XMM-Newton 
photometry, about 24% of the ASCA counts of source 19 could be due to confusion. We therefore 
excluded the ASCA data of 1993 and 2000 from the lightcurve shown in Figure SI On the other 
hand, it was noted that the luminosity during the high state in 1993 as suggested by the ASCA 
data, after subtraction of the flux contribution from confusion, at Lx ~ 2.9 x 10 39 erg s -1 would be 
similar in magnitude to the high-level Chandra flux in 2006 with Lx ~ 3. 1 x 10 39 erg s -1 . 

Based on the results above, we proposed that source 19 had undergone multiple state transi- 
tions. Starting from a low state in the 1991 ROSAT observation, it changed to the high/soft state in 
the ASCA observation in 1993, changed again to the VH state in the follow-up ASCA observation 
in 2000 and remained in this state until the first XMM-Newton observation in 2001 February. It 
then returned to the high/soft state in 2004 February and changed to the low/hard state six months 
later in 2004 August. In 2005 February, its luminosity was more than double the intensity level of 
1993 and stayed at the high/soft state through to the 2006 Chandra observation. If source 19 was 
truly in the VH state during the ASCA observation in 2000 and the XMM-Newton o bservation in 
2001, it would be the longest period that a VH state was observed (IBauer et al.ll2003|) . Finally, we 
searched for short-term variability within each XMM-Newton data but found none. 

The spectral/intensity variability of source 19 suggests this source is most likely a compact 
accreting object of a black- hole binary. A plausible scenario for the ULX source 19 is that it is 
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an IMBH formed via mergers of ma ssive stars/BH in a compact star cluster, as supported by the 
Comptonized disk model spectrum (Kubota et al.l 120021) . On the other hand, we cannot rule out 
other possible explanations for the su per-Eddington luminosity, such as a s tellar mass black hole 
with strongly beamed X-ray emission ( King et ahlboOl : Kording et al.120021) . Moreover, source 19 
was associated with a "tooth "-shaped optical nebula (IPakull & Mirionil2002|) . lRoberts et al.1 (12003b 
found two regions of O III emission located on the inside of the shell of the nebula and suggested 
they were caused by photonizat ion of the nebula shock excited from the ULX source 19. This 
scenario was later confirmed by lAbolmasov et al.l (|2007r ). LFeng & Kaaretl (|2008l) suggested that the 
incomplete shell in optical morphology could be a jet emission if the nebula was powered by an 
unusual powerful explosion in which a black hole was formed in source 19. 



5.2.2. Source 25 (IC342 X-2): An ultraluminous compact X-ray source 



This source is the second brightest ULX in IC 342 discovered by the Emsfem tlFabbiano & Trinchieri 



19871) . The spectrum taken by ASCA in 1993 show ed that it was in the low state, with absorbed 
^0.5-10 keV ~ 1.6 x 10 39 ergs _1 scaled to 1.8 Mpc (|Kubota et al.ll2001l ). and could be equally sat- 
isfactorily fitted by both power-law and disk blackbody model. The best- fit absorbed p ower-law 
parameters were N H = (14.3 ± 1.6) x 10 21 cm -2 and V = 1.39 ±0.10 (fokada et all 19981) . The sec- 
ond ASCA observation in 2000, with absorbed Lo.5-10 keV ~ 2.8 x 10 39 erg s -1 scaled to 1.8 Mpc, 
revealed that the spectrum of source 25 was more convex than in 1993 and could be expressed 
with a disk blackbody model of N n = (18.2±0.8) x 10 21 cm~ 2 and 7| n = 1 .62 ± 0.04 keV, whereas 
the power-law fit was unacceptable (iKubota et al.ll200ll) . The flux increased by a factor of ~ 2 
between these two observations. Such a transition from low/hard state to high/soft supported the 
black hole interpretation of source 25. The four XMM-Newton spectra of source 25 could be satis- 
factorily fitted by both the power-law and disk blackbody (diskbb) models, with the diskbb model 
giving consistently better fits than power-law models (Figure [6[ top-right). This was consistent 
with the spectral fits of the 1993 ASCA data for the source. The spectra had soft excess below 1 
keV and was the flattest in 2004 August with T = 1.3. In addition, the spectrum in 2004 August 
was different from the other three in which the inner disk temperature (7j n ) from diskbb model was 
particularly high together with a low photon index. Therefore, we also attempted the combined 
power-law and diskbb model and found that it gave a better fit to the soft excess seen below 1 keV 
with 7i n = 1-33 ±0.14 keV and T = 1.15 ±0.18. Similarly, a combined model of power-law and 
mekal gave similar satisfactory fit with kT = 0.16 ±0.05 keV and T = 1.54 ±0.05. Both models 
suggest the source was soft in 2004 August. These observations indicated that the source was 
undergoing significant spectral change from 2001 to 2005. 



Source 25 is the most variable source in our XMM-Newton catalog with a variability fac- 
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tor Sfl ux « 60. It became very luminous in the 2004 February data (unabsorbed Lo.5-10 keV ~ 
4 x 10 39 erg s -1 , even brighter than source 19) and in the 2006 Chandra observation (unabsorbed 
^0.5-10 keV ~ 3 x 10 39 erg s -1 ), implying an increase of a factor of > 4 in flux from its faintest state. 
With the low state being caught (2004 August and 2005 February) in between the recorded high 
states, the timescale of the high/soft-low/hard transition is much shorter for source 25 than for 
source 19. The intrinsic spectral/intensity change of source 25 is high enough to explain the transi- 
tion between the 1993 and 2000 ASCA observations and thus the contamination of nearby sources 
did not notably affect the observed variability if comparable magnitude of spectral/intensity tran- 
sitions also occurred in the ASCA data. Besides long-term variability, source 25 also exhibited 
short-term variability wit h a periodicity of 3 1 hr or 41 hr found in a long (~ 250 ksec net exposure 
time) ASCA observation (ISugiho et al.ll200lb . 



Similar to source 19, source 25 also suffered from the confusion problem with a large extrac- 
tion radius of 3' being used. In our XMM-Newton images, sources 23 and 34 in the ASCA analysis 
were included in the r = 3' circle centered on source 25, implyin g that the flux measured in the 
ASCA observations could be over-estimated by ~ 10%. However. iKongl (|2003|) noted that during 
the high/soft state in 2000, the source 25 was asymmetrically extending towards the direction of 
source 36, which was identified as a SSS located ~ 3' away. This suggested that the flux estimates 
of source 25 could also be affected by source 36. Source 36 is highly variable with flux changing 
by a factor of 5 between the observed faintest and brightest state and its brightness was similar to 
that of source 25 in 2004 August. In addition, the asymmetry of source 25 in 2000 might imply 
that source 36 was in its bright s tate, thus co ntributing more soft photons to result in a softer and 
brighter spectrum in source 25 (IKongl |2003|) as well as the spectral variability observed during 
1993 and 2000. Assuming nearby sources contributed also to the flux of source 25 in the ASCA 
observations in 1993 and 2000, the source could be dimmed by more than a factor of 2 from its 
high state in 2000 (with unabsorbed Lo.5-10 keV ~ 2 x 10 39 erg s -1 ) to low state in 2001 (with un- 
absorbed Lo.5-10 keV ~ 1 x 10 39 erg s -1 ). All these supported the proposition that source 25 is an 
accreting binary object. Nevertheless, there is a lack of multiwavelength analysis for this source. 
Studies of optical or radio counterparts will help to identify the nature of the compact object. 



5.2.3. Source 38 (IC342 X-3): The nuclear X-ray source 

Source 38 is one of the three historical ULXs detected by Einstein. Recent high resolution 
X-ray imaging studies show ed that the sour ce was associated with the galactic center and thus 
confirmed it was not a ULX (Ma k et al A comprehensive study of the multiwavelength and 



spatial analysis of this source based on the Chandra HRC-I observation in 2006 was also presented 
in that paper. Here we focus on the spectral analysis based on the four XMM-Newton observations. 
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We found that simple models su ch as power- l aw, M EKAL, Raymond-Smith, and MCD all 
gave unacceptable fits, with \v > 1.7. Bauer et al.l (|2003r) obtained a good fit for source 38 using 



the XMM-Newton data in 2001 with a best fit model of an absorbed MEKAL + power-law model 
(kT = 0.30^$ keV, V = 2.52+g;}|, N H = 6.4+^). We therefore attempted to fit source 38 with 
this model and with a combination of different absorption models. We found a combination of 
MEKAL and power-law model with Galactic absorption gave the best fit. Parameters of the spectra 
fit s were listed in Tab le M The best fit A^h was about (3.5-7.2) x 10 21 cm~ 2 , similar to the result 



of iBauer et al.l (|2003f ). The best-fit photon index T showed a general decreasing trend from 2001 
February (2.6) to 2005 February (1.5). Except for the 2004 February data (V = 2.2), the spectra 
of source 38 were outside the typical range for AGNs (T = 1.7-2.3). In addition, there were 
also trends of decreasing absorption column density, and increasing plasma temperature (kT) from 
2001 to 2005, strongly supporting a true spectral change. 

Closer inspections of the spectral shapes over the 4 years of XMM-Newton data also revealed 
that the spectrum in 2001 was very different from the others (Figure [6] bottom-right). There was 
an abrupt change in the slope for emission above 1 keV which showed a much harder emission in 
2001 than in the follow-up observations. On the other hand, emission lines at ~0.8 keV (Fe L), 
1.1 keV (Ne Lyr), 1.4 keV (Mg XII), and 1.9 keV (Si XIII) were most prominent in the spectra of 
both 2004 observations, primarily because of the higher signal-to-no ise of these data sets. These 
emission lines were also observed in the nucleus of NGC 1808 (see Jimenez-Bailon et al. 2005i 



fig. 6), and we adopted their line identifications here. No Fe Kr line was detected with high 
significance. This motivated us to replace the MEKAL model with one in which the abundances 
of individual elements could be fitted, that is, the VMEKAL model. The abundances of O, Ne, 
Mg, Si, and Fe were left as independent free parameters in the fits while the abundances of other 
elements are set to solar values. The best-fit abundances were sub- solar, except for the 2004 
February observation in which an abundance of > 18 Z , much higher than the value obtained 
by B03 (Z = 6.39Z Q for MEKAL component) was found. The values of the best-fit x 2 /dof of 
the VMEKAL and MEKAL models were compared independently in each observation using F- 
test. The F-statistics suggested that the \ 2 V improvement obtained when replacing MEKAL with 
VMEKAL was statistically significant only for the two highest S/N 2004 observations. In addition, 
an intrinsic absorption to the power-law component did not improve the fit and this intrinsic column 
density Nu gave a value of zero. 

An alternative model, with a power-law + blackbody wa s also compa tible with the data, and 



had been used to fit the spectrum for the 2001 observation by iKongl (|2003|) (kT = 0. 1 1 keV, A^h = 
8.7 x 10 21 cm~ 2 , T = 2). Even after we revised the model by adding a gaussian line at ~ 0.8 keV to 
provide best fits there were large deviation in the spectra fits through the whole spectrum, and the 
power-law + MEKAL offered better fits to the data. Still, this model gave a reasonable estimate to 
the luminosity of the source. 
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A final remark is that the variation of the spectral fits of source 38 to both power-law and 
MEKAL models is very consistent. It is interesting to note that spectral change seen in nuclear X- 
ray source is very rare and the change seen in source 38 would be supplement to this rare sample. 
On the other hand, the underlying mechanism that drives this spectral change in galactic nuclear 
source is still unclear. Further investigations are needed to understand the emission mechanism. 



5.2.4. Source 27: A Supernova Remnant candidate 



As described in Section l4~3.2[ source 27 had been identified as a supernova remnant associated 
with the SNR object 1 of IPodorico et al.lll980l We could not identify any X-ray sources from 
previous observations with Einstein, ASCA, and ROSAT. It was first detected in the 2001 February 
XMM-Newton observation and was detectabl e in all observati ons at Lx ~ 1 x 10 37 erg s -1 until the 
last observation in 2006 by Chandra (C7 in iMak et al.ll2008|) . While the detection sensitivity of 
Einstein and ASCA were an order of magnitude worse than XMM-Newton and Chandra, ROSAT 
had comparable sensitivity as these two satellites and had observed the field of source 27 in 1991. 
We used PEVIMS to estimate the ROSAT PSPC count rates at the source position. Assuming 
the source flux had not changed significantly over the years, and using a power-law model with 
T = 2 and Nu = 8 x 10 21 cm -2 , the deduced ROSAT count rate for source 27 would have been 
3.8 x 10~ 4 cps in the 0. 1 -4 keV band. We compared this predicted count rate with the background 
subtracted source counts of the ROSAT data at the source 27 position using DMEXTRACT. The 
measured ROSAT count rate (1 .9 x 10~ 4 cps) is found to be at least a factor of 2 below the expected. 
This could suggest a possib le scenario of an X-ray bin ary associated with the SNR that turns the 
source off and on again (e.g JWilliams et al.ll2005L l20Q7f) . This could be confirmed with a long-term 
monitoring of the source in the future. This SNR was described as a diffuse shell with an angular 
size of 42 pc (assumed distance of 2.9 Mpc) in the optical images (|Dodorico et al.M 198QQ . After 
scaling to our assumed distance of 1.8 Mpc, this corresponds to a size of ^16 pc. However, source 
27 was not resolved in the high spatial resolution Chandra HRC-I data and was consistent with a 
point source when compared to the Chandra PSF of size 9 w 2" (corresponding to 18 pc). This 
X-ray spatial extent is comparable to that in the optical images. It is worth noting that there are 
only a few sp atially resolved X-ray SNRs beyond the Milky Way and Magellanic Clouds observed 
recently (e.g. Kong et alj|2002al : iKong et al.ll2003l : iKong et alj|2004bl) . and the unresolved X-ray 
structure of source 27 is possibly due to its intrinsic size being smaller than the Chandra PSF. 

We noted that the Raymond-Smith (RS) and nonequilibrium ionization (NEI) models were of- 



ten a pplied to study X-ray emission from extragalactic SNRs (e.g. IKong et all2002al : ISchlegel et al 



2000|) . The RS model is a simple collisional equilibrium ionization model, while the NEI model is 
appropriate for modeling SNRs whose ages are smaller than the time required to reach ionization 
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equilibrium. The NEI model consists of an electron temperature (kT e ) and an ionization timescale 
(n s t), where n e and t are the mean electron density and the elapsed time after the plasma is shock 
heated to a constant temperature kT s respectively. Due to limited counts, we were only able to fit 
the spectra of the two 2004 observations. The NEI model gave much better fits to the data than 
other simple models (including RS) since these simple models gave only acceptable fits to the 
spectrum (y} v > 2) in 2004 August. Results of the spectral fits, assuming solar abundance, were 
given in Tableland Figure [6] (bottom- left). We attempted to fit the abundances for the NEI model 
but could not constrain the parameters. Moreover, there are indications of line emissions at 0.9 and 
2 keV. It is likely that the 0.9 keV feature comes from the Fe L shell lines and Ne K shell lines. The 
fitted electron temperature at above 2 keV was relatively high but it had been also observed in some 
extragalactic SNRs (e.g., N132D in LMC LFavata et al.lll997|) which could be d ue to a shock-heate d 
swept-up circumstellar mediu m or a inhomogeneity of the interstellar medium (|Kong et al ] l2004bh . 
Following iKong et al.l (|2002ar ). we estimated the physical parameter s of the SNR th rough the Se- 
dov solution, assuming an initial explosion energy of 3 x 10 50 ergs (|Blair et al. 1981 ). a radius of 9 
pc, and an electron temperature of 2.37 ± 1.24 keV (average value from the two spectral fits). We 
derived the age of the SNR to be 2500^qq° years and the number density of the ambient gas to be 
0.10^[| q3 cm -3 . The derived age and number density are both relatively small compared to other 
X-ray observed SNRs. This could have explained why this SNR was not detected in the X-ray in 
previous missi ons as nondetected X -ray SNRs usually reside in regions with ambient densities less 
than 0.1 cm" 3 dlVIagnier et all 1997b . 
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Table 1. 



Archival XMM-Newton Observations of IC342 







Aim point 




T exp (ksec) 




Obs ID 


Date 


RA (J2000.0) 


DEC (J2000.0) 


EPIC PN 


EPIC MOS1 


EPIC MOS2 


0093640901 


2001-02-11 


03:46:49.40 


+68:05:38.9 


5.6 


9.5 


9.6 


0206890101 


2004-02-20 


03:46:15.92 


+68:08:43.6 


13.0 


20.1 


20.7 


0206890201 


2004-08-17 


03:45:56.41 


+68:07:23.7 


19.7 


23.6 


23.5 


0206890401 


2005-02-10 


03:45:55.66 


+68:07:31.8 


6.1 


6.1 


5.9 



Table 2. XMM-Newton X-ray source catalog of IC342 



ID 


Source Name 


RA 


DEC 


Positional error 




Net counts 




HR1 


HR2 


L x (0.3-12 keV) 


Note 




XMMU J 


(J2000.0) 


(J2000.0) 


(arcsec) 


soft (0.3-1 keV) 


medium (1.0-2.0 keV) 


hard (2.0-12.0 keV) 






(xlO 38 ergs" 1 ) 




1* 


034333.4+680224 


03:43:33.47 


68:02:24.8 


1.34 


88.9 ±10.0 


181.2± 14.0 


271.3± 19.5 


0.34 ±0.07 


0.51 ±0.07 


0.12 ±0.02 




2* 


034354.8+675657 


03:43:54.88 


67:56:57.4 


1.41 


733.1 ±31.9 


292.4 ±19.8 


83.3 ±29.5 


-0.43 ±0.04 


-0.80 ±0.07 


0.67 ± 0.06 




3 


034402.1+680935 


03:44:02.13 


68:09:35.8 


0.73 


124.4 ±13.3 


391.7±21.8 


434.5 ±27.7 


0.52 ±0.06 


0.55 ±0.06 


0.55 ±0.04 




4 


034426.8+681739 


03:44:26.82 


68:17:39.0 


0.85 


195.0±15.2 


74.9±9.8 


< 8.9 


-0.45 ±0.07 


-0.75 ±0.11 


0.77 ±0.05 


f ,t 


5 


034441.1+680827 


03:44:41.10 


68:08:27.7 


0.83 


< 21.8 


114.9 ±13.0 


166.0±18.8 


0.95 ±0.17 


0.97 ±0.16 


0.10±0.01 




6 


034447.4+680840 


03:44:47.40 


68:08:40.8 


0.63 


85.2± 12.1 


426.6 ±22.3 


503.2 ±27.2 


0.67 ±0.06 


0.71 ±0.06 


0.50 ±0.05 


f , 


7 


034449.0+680736 


03:44:49.06 


68:07:36.3 


1.34 


74.4 ±9.2 


< 17.2 


< 17.2 


-0.72 ±0.15 


-0.82 ±0.17 


0.31 ±0.04 


f ,t 


8 


034449.6+680216 


03:44:49.65 


68:02:16.3 


0.80 


23.7±6.8 


80.6 ±10.6 


105.0±14.3 


0.55±0.14 


0.63 ±0.15 


0.16 ±0.03 




9 


034504.5+675910 


03:45:04.56 


67:59:10.7 


1.49 


<26.5 


18.8±5.7 


33.4±8.8 


0.38 ±0.27 


0.60 ±0.27 


0.09 ±0.03 


f , s 


10 


O345O7.5+680111 


03:45:07.51 


68:01:11.7 


0.60 


102.0±10.7 


169.5 ±13.8 


115.3± 12.7 


0.25 ±0.07 


0.06 ±0.08 


0.74 ±0.07 


f ,t 


11 


034509.3+680855 


03:45:09.32 


68:08:55.7 


0.80 


32.7 ±8.6 


84.5 ±11.3 


80.0 ±16.0 


0.44 ±0.1 3 


0.42 ±0.17 


0.09 ±0.02 


f , 


12 


034510.1+680231 


03:45:10.12 


68:02:31.0 


0.44 


253.5 ±17.1 


956.5±31.8 


875.8 ±32.4 


0.58 ±0.03 


0.55 ±0.04 


0.87 ±0.05 


f ,t 


13 


034515.3+681648 


03:45:15.37 


68:16:48.5 


0.95 


< 11.5 


23.7±8.2 


93.2±13.9 


0.44 ±0.32 


0.82±0.19 


0.21 ±0.03 


f,t 


14 


034531.3+681912 


03:45:31.32 


68:19:12.0 


0.68 


123.3 ±13.9 


520.0 ±24.8 


584.1 ±31.2 


0.62 ±0.05 


0.65 ±0.06 


0.72 ±0.05 




15 


034533.9+675912 


03:45:33.93 


67:59:12.9 


0.51 


63.9±8.9 


304.7 ±18.6 


489.0 ±24.5 


0.65 ±0.07 


0.77 ±0.06 


0.36 ±0.03 




16 


034536.9+680733 


03:45:36.93 


68:07:33.4 


0.90 


15.3 ±6.1 


53.0±9.1 


37.0 ±10.7 


0.55 ±0.18 


0.42 ±0.25 


0.07 ±0.02 


s 


17 


034539.9+680309 


03:45:39.91 


68:03:09.4 


0.29 


330.2 ±19.4 


1230.0 ±36.2 


1789.6 ±44.6 


0.58 ±0.03 


0.69 ±0.03 


1.56 ±0.07 


f 


18 


034547.1+680547 


03:45:47.12 


68:05:47.5 


0.61 


< 34.3 


51.2±10.4 


71.6± 13.6 


0.74 ±0.25 


0.81 ±0.24 


0.10±0.01 


f ,t 


19 


034555.6+680455 


03:45:55.62 


68:04:55.9 


0.09 


3846.8 ±65.4 


26206.4 ±163.9 


32480.2 ±183.9 


0.74 ±0.01 


0.79 ±0.01 


21.00±0.18 


f,s 


20 


034556.8+675925 


03:45:56.86 


67:59:25.8 


0.76 


45.6±7.9 


91.3 ±11.1 


156.5 ±15.5 


0.33 ±0.11 


0.55 ±0.10 


0.19 ±0.03 


f,t,s 


21 


034559.4+680537 


03:45:59.44 


68:05:37.6 


0.29 


105.9 ±12.3 


555.2±26.7 


799.6 ±3 1.9 


0.68 ±0.05 


0.77 ±0.05 


0.83 ±0.04 




22 


034602.0+675513 


03:46:02.09 


67:55:13.7 


1.11 


< 30.7 


20.9 ±6.1 


27.5 ±8.9 


0.35 ±0.26 


0.47 ±0.29 


0.14 ±0.02 




23 


034606.1+681029 


03:46:06.19 


68:10:29.3 


0.43 


41.0±8.6 


228.6 ±16.8 


286.4 ±22.1 


0.70 ±0.09 


0.75 ±0.09 


0.30 ±0.02 


s 


24 


034606.5+680705 


03:46:06.54 


68:07:05.3 


0.33 


976.2 ±32.1 


358.8 ±19.8 


48.3 ±10.5 


-0.46 ±0.03 


-0.91 ±0.04 


3.24±0.10 


f,s 


25 


034615.6+681112 


03:46:15.64 


68:11:12.2 


0.12 


513.4±28.7 


9198.3 ±98.2 


41562.4 ±207.3 


0.89 ±0.01 


0.98 ±0.01 


12.0()±0.11 


f,s 


26 


034619.5+680554 


03:46:19.59 


68:05:54.7 


0.83 


26.0 ±6.0 


25.4±6.8 


12.9±7.8 


-0.01 ±0.18 


-0.34 ±0.27 


0.10 ±0.02 


f 


27 


034622.1+680506 


03:46:22.12 


68:05:06.9 


0.55 


42.2±7.9 


262.5 ±17.0 


154.6±14.3 


0.72 ±0.08 


0.57 ±0.10 


0.18 ±0.02 


s 


28 


034625.9+680421 


03:46:25.91 


68:04:21.8 


0.75 


30.1 ±6.1 


24.1 ±5.7 


<22.0 


-0.11 ±0.16 


-0.43 ±0.21 


0.10 ±0.03 


f 


29 


034626.7+680455 


03:46:26.70 


68:04:55.4 


0.75 


52.6±7.9 


42.1 ±7.5 


18.1 ±6.6 


-0.11 ±0.12 


-0.49 ±0.16 


0.18 ±0.03 


f,t 


30 


034627.5+680410 


03:46:27.56 


68:04:10.7 


0.82 


27.8±5.7 


37.2±6.5 


<7.2 


0.14±0.13 


-0.61 ±0.24 


0.19 ±0.03 


f,t,s 


31 


034633.4+675855 


03:46:33.48 


67:58:55.8 


0.87 


31.7±7.0 


98.9 ±10.9 


139.1 ±14.7 


0.52 ±0.11 


0.63 ±0.11 


0.11 ±0.02 




32 


034638.9+675551 


03:46:38.94 


67:55:51.6 


0.79 


53.9±8.4 


158.7±13.4 


106.0±13.0 


0.49 ±0.08 


0.33 ±0.10 


0.38±0.01 


f,s 


33 


034639.8+680517 


03:46:39.80 


68:05:17.5 


0.59 


<9.6 


<7.9 


256.5 ±19.0 


0.00 ±0.00 


1.00 ±0.11 


0.23 ±0.03 


f 


34 


034643.6+680611 


03:46:43.65 


68:06:11.4 


0.33 


< 12.9 


21.3±7.5 


49.1 ±10.3 


0.49 ±0.37 


0.74 ±0.26 


1.22 ±0.05 


s 



Table 2 — Continued 



ID 


Source Name 

~V~\ f~\ AT T T 

XMMU J 


RA 
(J 2000.0) 


DEC 
(J 2000.0) 


Positional error 
(arcsec) 


sort (0.3-1 kev) 


Net counts 
medium (1.0-2.0 kev) 


hard (2.0-12.0 kev) 


HR1 


HR2 


Lx (0.3-12 kev) 

i — i n38 ^—1 \ 

(x 10 ergs ) 


Note 


35 


034644.0+681104 


03:46:44.07 


68:11:04.6 


0.80 


189.3±22.1 


729.3 ±37.2 


812.1 ±38.9 


0.59 ±0.05 


0.62 ±0.05 


0.07 ±0.01 


f 


36 


034645.1+680946 


03:46:45.14 


68:09:46.1 


0.22 


2045. 2 ±46.1 


1642.4 ±41. 5 


307.2 ± 18.9 


-0.11 ±0.02 


-0.74 ±0.03 


4.93 ±0.10 


f , t , s 


37 


034646.1+680524 


03:46:46.19 


68:05:24.7 


0.22 


62.2±8.6 


63.7 ±9.1 


< 9.5 


0.01 ±0.10 


-0.96 ±0.22 


1.22 ±0.07 


f , t 


38 


034648.5+680547 


03:46:48.52 


68:05:47.4 


0.20 


3517.3 ±63.0 


393 1.4 ±68.0 


1497.1 ±45.4 


0.06 ±0.01 


-0.40 ±0.02 


9.63±0.16 


s 


39 


034651.3+680028 


03:46:51.32 


68:00:28.6 


0.54 


36.4±8.3 


277. 6± 18.2 


477.2 ±24.8 


0.77 ±0.08 


0.86 ±0.07 


0.37 ±0.03 


f , 


40 


034652.5+680535 


03:46:52.59 


68:05:35.4 


0.46 


64.1 ± 12.3 


74.6 ± 13.9 


49.3 ± 13.4 


0.08±0.13 


-0.13 ±0.16 


1.07 ±0.07 


f , t , s 


41 


034652.6+680847 


03:46:52.63 


68:08:47.2 


0.64 


< 21.2 


67.9±9.7 


130. 0± 13.8 


0.86 ±0.20 


0.92 ±0.15 


0.10±0.01 


f , t 


42 


034652 7+680504 


03-46-52 74 


68-05-04 9 

UO.UJ An. y 


0.69 


43 9-1-9 9 


1 <0 4 + 1 4 7 


213.4-1- 17.4 


"j's + n 10 


66 ± 09 


?6 ± 03 


f t 


43 


034654.3+675901 


03:46:54.36 


67:59:01.2 


0.86 


38.7±7.2 


80.5 ±10.0 


<23.7 


0.35 ±0.11 


-0.56 ±0.23 


0.17 ±0.03 


f,t,s 


44 


034657.2+680619 


03:46:57.27 


68:06:19.4 


0.23 


495.6 ±23.8 


1807.4 ±43.9 


1909.0±45.4 


0.57 ±0.02 


0.59 ±0.02 


2.93 ±0.08 


f 


45 


034659.3+680316 


03:46:59.31 


68:03:16.0 


0.54 


300.1 ±18.4 


914.7±31.4 


565.7 ±26.0 


0.51 ±0.03 


0.31 ±0.04 


0.80 ±0.03 


f,t,s 


46 


034701.9+680237 


03:47:01.93 


68:02:37.6 


0.52 


83.7± 11.9 


261.0±18.3 


359.2±21.9 


0.51 ±0.07 


0.62 ±0.07 


0.44 ±0.03 




47 


034703.8+680905 


03:47:03.83 


68:09:05.0 


0.58 


68.5 ±10.0 


151.7±13.3 


239.9 ±17.4 


0.38 ±0.08 


0.56 ±0.07 


0.30 ±0.02 




48 


034704.6+680517 


03:47:04.69 


68:05:17.4 


0.68 


141.6±13.3 


102.4±12.0 


<35.7 


-0.16 ±0.07 


-0.70 ±0.12 


0.33 ±0.03 




49 


034708.1+681059 


03:47:08.12 


68:10:59.0 


0.64 


88.0±10.1 


299.2±18.2 


231.5 ±17.7 


0.55 ±0.06 


0.45 ±0.07 


0.69 ±0.06 


f,t 


50 


034718.6+681128 


03:47:18.63 


68:11:28.5 


0.51 


151.4±13.8 


455.4 ±22.6 


468.9 ±25.0 


0.50 ±0.05 


0.51 ±0.05 


0.79 ±0.05 


f 


51 


034722.3+681538 


03:47:22.38 


68:15:38.5 


0.84 


61.0±8.5 


46.6±7.9 


< 8.6 


-0.13±0.11 


-0.64 ±0.1 8 


0.31 ±0.04 


f,t,s 


52 


034722.9+680859 


03:47:22.96 


68:08:59.5 


0.40 


317.2±19.0 


1033.8 ±33.1 


375.5 ±22.7 


0.53 ±0.03 


0.08 ±0.04 


0.95 ±0.05 


f,t,s 


53* 


034726.5+680849 


03:47:26.58 


68:08:49.2 


0.94 


<() 


<() 


88.9 ±10.4 


0.07 ±0.30 


0.81 ±0.15 


0.12 ±0.02 


f,t 


54 


034748.1+681530 


03:47:48.12 


68:15:30.1 


0.91 


23.4±6.8 


49.2±8.8 


41.3±12.8 


0.35 ±0.16 


0.28 ±0.23 


0.11 ±0.03 


f,t 


55 


034803.5+681114 


03:48:03.58 


68:11:14.9 


1.20 


<6.4 


36.1 ±7.9 


57.0±11.7 


0.93 ±0.32 


0.95 ±0.29 


0.04 ±0.02 


s 


56 


034805.1+680137 


03:48:05.10 


68:01:37.9 


1.37 


27.6±6.1 


16.9±4.9 


< 8.1 


-0.24 ±0.1 8 


-1.00±0.41 


0.15 ±0.04 


f,t,s 


57 


034807.1+680455 


03:48:07.10 


68:04:55.5 


0.34 


329.3 ±22.2 


1202.9 ±39.4 


1857.1 ±52.0 


0.57 ±0.03 


0.70 ±0.03 


1.57 ±0.07 


f , 


58 


034817.9+680204 


03:48:17.99 


68:02:04.4 


0.99 


44.9 ±11.1 


191.9± 19.3 


399.1 ±30.6 


0.62 ±0.11 


0.80 ±0.09 


0.36 ±0.02 





Table 2 — Continued 



ID 


Source Name 


RA 


DEC 


Positional error 




Net counts 




HR1 


HR2 


L x (0.3-12 keV) 


Note 




XMMU J 


(J2000.0) 


(J2000.0) 


(arcsec) 


soft (0.3-1 keV) 


medium (1.0-2.0 keV) 


hard (2.0-12.0 keV) 






(xlO 38 ergs" 1 ) 




59 


034825.4+680817 


03:48:25.42 


68:08:17.1 


1.46 


19.4±5.5 


38.3±9.0 


85.6±14.3 


0.33±0.19 


0.63±0.17 


0.14 ±0.03 


f ,t 


60 


034832.3+675644 


03:48:32.37 


67:56:44.3 


2.00 


< 9.3 


17.1 ±6.6 


64.0 ±10.6 


1.00 ±0.63 


1.00±0.25 


0.23 ±0.02 




61 


034914.5+675732 


03:49:14.58 


67:57:32.4 


1.23 


19.5±7.0 


49.8 ±11.7 


242.4±17.5 


0.44 ±0.21 


0.85 ±0.09 


0.19 ±0.05 





Note. — Column (1) lists the running ID number of the X-ray sources. Column (2) lists the X-ray source names given as XMMU JHHMMSS.S+DDMMSS. Column (3) and (4) list the X-ray source positions 
that used the astrometric reference as discussed in Section [3~T1 Column (5) lists the rms positional errors returned by the WAVDETECT task after aligning the four observations. The cross registration errors of 
aligning the observations and the astrometric registration errors with the USNO B1.0 catalog are not included. Column (6)-(8) list the total background and exposure corrected net counts and its errors of X-ray 
sources for all 4 observations in the soft, medium and hard band respectively, as determined by DMEXTRACT. For sources that have S/N < 2.5 in a particular band, the 3a upper limit are used instead. Column 
(9) and (10) list the average HR1 and HR2 values, defined as HR1 = (M — S)/(M + S) and HR2 = (H — S)/(H + S), and their corresponding errors averaged from data sets where the the source was present. 
Column (11) lists the unabsorbed 0.3 — 12 keV luminosities in unit of (x 10 38 erg s ) determined by EMLDETECT, by assuming an absorbed power-law spectrum with r = 2 and Nu = 8 x 10 21 cirT 2 , averaged 
from data sets where the source was present. A distance of 1.8 Mpc to IC342 is assumed. Column (12): Symbols are defined as follows: t = X-ray transient; f = flux variable; s = spectral variable. 
* These sources were only detected with EDETECT_CHAIN software of SASDAS and not by the WAVDETECT task of CIAO. 



Table 3. Cross -correlation of our XMM-Newton catalog with Chandra, ROSAT, Einstein, and 

previous XMM-Newton catalogs 



Source ID 


previous XMM-Newton " (offset ") 


Chandra b (offset ") 


ROSAT c (offset ") 


Einstein d (offset ") 


6 


XI (3.2") 


- 


- 


- 


7 
8 


X3 (5.6") 
X2 (1.6") 


- 
- 


- 
- 


- 
- 


9 


X4 (4.2") 


- 


- 


- 


12 


- 


CI (5.7") 


- 


- 


15 


X5 (3.2") 


- 


- 


- 


17 


X6(1.8") 


C2 (3.5") 


Rl (1.0") 


- 


19 


X7 (2.3") 


C3 (0.1") 


R3 (1.6") 


IC342 X-l (25.2") 


20 


X8 (2.6") 


- 


- 


- 


21 


X9 (3.4") 


C4(l.l") 


- 


- 


22 


X10 (4.2") 


- 


- 


- 


23 


XI 1 (2.3") 


- 


- 


- 


24 


X12(1.2") 


C5 (1.2") 


R4 (0.4") 


- 


25 


X13 (3.7") 


C6 (2.4") 


R5 (7.3") 


IC342 X-2 (30.6") 


27 


X14(3.1") 


C7 (0.9") 


- 


- 


28 


X15 (6.8") 


- 


- 


- 


29 


X16 (2.4") 


- 


- 


- 


30 


- 


C8 (3.2") 


- 


- 


31 


X17 (3.6") 


- 


- 


- 


32 


X18 (3.7") 


- 


- 


- 


35 


X19 (2.0") 


C9 (0.7") 


R6(1.6") 


- 


36 


X20 (3.5") 


Cll (2.0") 


R7 (5.2") 


- 


38 


X21 (1.7") 


C12(l.l") 


R8 (3.1") 


IC342 X-3 (40.7") 


jy 


A.ZZ \L.L ) 








41 


X23 (4.8") 








42 


X24 (2.1") 


C14(l.l") 






43 


X25 (5.0") 








44 


X26 (2.1") 


C15 (1.1") 


R9 (4.5") 




46 


X27 (0.9") 


C16(0.7") 






47 


X28 (1.5") 


C17 (3.2") 






48 


X29 (3.5") 


C18(1.2") 






50 


X30 (3.4") 


C21 (0.8") 






52 


X31 (2.2") 


C22(l.l") 






56 


X32 (5.7") 








57 


X33 (3.1") 


C23 (2.9") 


R10 (6.5") 




58 


X34 (7.1") 








59 


X35 (3.9") 








60 


X36 (5.5") 








61 


X37 (0.9") 









Note. — " T able 1 of lKond |2003|) .I Bauer et al ] 120031) also gave the analysis of the same XMM-Newton observation 
as lKondl2003l). but the catalo gs are similar and lKond |2003|) has two more sources. Searching radius = 6" 

Table 1 o f Mak et al.| |20 08jl. Pre fix "C" denotes Chandra sources. Searching radius = 6" 
c Table 1 o flB^grriarietaiTI 19931). Prefix "R" denotes ROSAT sources. Searching radius = 6" 

d Table 3 of Fabbiano & Trinchieril h987l) . Prefix "IC342-X" denotes the 3 brightest ULX in IC342 identified by Einstein. 
Searching radius = 40" 



Table 4. Energy Conversion Factors for each instrument and energy band 



Detector 


Soft (0.3-1 keV) 


ECF(10- U ) 
Medium (1-2 keV) 


Hard (2-12 keV) 


MOS1 


0.032 


0.597 


0.266 


MOS2 


0.032 


0.597 


0.266 


PN 


0.0987 


1.24 


0.516 



Note. — Count rate to flux conversion factors for each EPIC instruments in 
the soft, medium, and hard band, assuming an absorbed power-law spectrum 
with a photon index of 2 and absorption column density of 8 x 10 21 cirT 2 . 



Table 5. Classification Scheme of the X-ray sources 



Object Type 


Hardness Ratio 


Catalog 


Selection criteria 
X-ray Spectral properties 


Others 


XRB 




a 


XRB: diskbb 1 with kT ss 1 -3 keV 


XRBs: variable 


sss 


c.f. fig. 3 in D03* 




bb with kT < lOOeV 




SNR 


HRl z > 0.1 and HR2 Z < -0.4 (Z06) c 


Do80 rf 


dominated by thermal spectrum < 2 keV 




Stars 


HR\ Z < 0.3 and HR2 Z < -0.4 


USNOBl.C 


RS 2 with kT m 1 keV 


J-K < 0.8 and J < 12.5, log(/ x // op t) < 1 


X-ray transient 








1. variability factor > 3 



2. at least detected in one observation 
with Lx > 10 37 erg s~' and not detected in 
at least one other observation 



Note. — 1 Disk Blackbody model; 2 Raymond-Smith model 

"XMM -Newton : lKonj2q03t Chandra: Mak, Pun. & Kond2008t ROSAT: iBregman, Cox. & Tomisakal 19931 ; Einstein: iFabbiano & Trinchierill 19871 
foi Stefano & Konj2003tt 

Misanovic et al. 2006; 

Dodorico et aljl98d ; 

Monet et alj2003l 



Table 6. Summary of the identified SSS (source 36) and QSS (sources 52, 54, 56) 



Source ID 


Observation 


Soft (0.1-1.1 keV) 


Medium (1.1-2 keV) 


Hard (2-7 keV) 


HR1 


HR2 


36 


2002 February * 


170.89 ±13.35 


6.81 ±2.87 


3.05 ±2.57 


-0.92 


-0.96 




2004 February 


627.06 ±25.44 


132.65 ±12.13 


18.82±5.77 


-0.65 


-0.94 




2004 August 


1738.74±42.56 


1047.46 ±33.09 


270.79 ±16.86 


-0.25 


-0.73 




2005 February* 


18.65±4.76 


1.48±2.13 


0.00 ±1.99 


-0.85 


-1.00 


52 


2002 February 


66.3 ±8.6 


125.6± 11.4 


51.1 ±7.7 


0.31 


-0.13 




2004 February 


216.1 ±15.6 


31 1.1 ± 18.5 


109.8 ±12.3 


0.18 


-0.33 




2004 August 


206.8 ±15.3 


434.8 ±21.4 


159.4± 13.5 


0.36 


-0.13 




2005 February * 


22.1 ±5.9 


14.2±4.5 


0.0±4.5 


-0.22 


-1.00 


54 


2002 February 


3.4±3.9 


0.0 ±2.9 


0.0±3.6 


-1.00 


-1.00 




2004 February* 


2.0±8.9 


42.0 ±9.6 


0.0 ±10.6 


0.91 


-1.00 




2004 August 


22.6 ±7.5 


48.8 ±8.6 


25.8±8.5 


0.37 


0.07 




2005 February* 


11.4±6.4 


6.6±5.1 


0.0 ±5.4 


-0.27 


-1.00 


56 


2002 February* 


26.8 ±6.3 


14.9 ±4.6 


0.0±3.4 


-0.29 


-1.00 




2004 February 


2.4±9.1 


2.4±7.6 


19.5 ± 11.0 


0.00 


0.78 




2004 August 


15.7±7.8 


10.8±7.2 


9.1±6.8 


-0.18 


-0.27 




2005 February 


3.3±5.4 


0.7 ±5.1 


1.2±7.4 


-0.64 


-0.47 



Note. — 



* Observations in which the source satisfied the selection criteria of SSS or QSS. 



Table 7. Identification of SSS, QSS, SNR, and foreground star 



Source ID 


Type 1 


Identification 2 


Radial Offset 3 
USNOB1.0 2MASS 


D08 2 


2 


Star 


Cat3 


4.0" 


3.9" 




7 


Star 


Catl 


4.0" 


0.7" 




16 


Star 


Cat2 


3.1" 


3.3" 




24 


Star 


Catl 


1.5" 


1.2" 




26 


Star 


Catl 


3.1" 


3.1" 




27 


SNR 


D80 






9.4" 


28 


Star 


Cat2 


4.3" 


5.2" 




29 


Star 


Catl 


1.7" 


2.2" 




36 


SSS 


HR 








37 


Star 


Cat2 


1.8" 






43 


Star 


Cat2 


2.5" 


2.2" 




48 


SNR 


HR 


0.6" 






51 


Star 


Cat2 


1.5" 


2.2" 




52' 


QSS-noh 


HR 








54* 


QSS-noh 


HR 








56 


QSS-noh 


HR 








61 


Star 


Cat2 


2.6" 


3.1" 





Note. — * These sources are classified as SSS/QSS candid ates only. 

1 Class ification of SSS and QSS adopted the scheme by Di St efano & Kond 
<2003bh . QSS-noh: sources exhibited little or no emission above 1 keV; 

2 Foreground star classification (Catl: USNOB1.0 + X-ray to optical flux ratio 
+ 2MASS + hardness ratio; Cat2: USNOB1.0 + X-ray to optical flux ratio + 
hardness ratio; Cat3 : USNOB1.0 + X-ray to optical flux ratio + 2Mass); D80: 
lDodoricoetalll98ri HR: hardness ratio criteria 

3 Searching radius of catalogs: USNOB1.0 = 5"; 2MASS = 5" D80 = 10" 



Table 8. Spectral fits to the brightest 1 1 X-ray sources in IC342 



Source ID 


Obs year 


XSPEC model 


N H (10 22 crrT 2 ) 


T/logKO" 


£7* (keV) 


Xl'dof 


Z>x (10 38 ergs-') 




Flux variables 1 


17 


2004 February 
2004 August 


pow 
diskbb 

pow 
diskbb 


0.20^ 
u - ;>;, -0.13 

A -3 A+0.08 

U - JU -0.07 


1 26 +03 ° 
1 - ZD -0.26 

i n q +0 2 ° 


7 70+I.I4 
z - /J -0.62 

1 7-r+0.26 


1.39/21 
1.16/21 
1.35/42 

1 1 S IA1 


0.73 
0.53 
0.73 

U.o / 


19* 


2001 February 

2004 February 
2004 August 

2005 February 


pow 
diskbb 

pow 
diskbb 

pow 
diskbb 

pow 
diskbb 

RS 


A <r 7 +0.06 
"■■"-0.05 

U ' J1 -0.03 

A QO+0.03 

U.S3_ 002 
U - 48 -0.01 

u - DZ -0.03 
35+O-02 
u ' JJ -0.02 

A Q 1+0.04 

U ' M -0.04 
a 4 7 +O.02 
U '^'-0.02 
a fi7 +0.03 
u -°'-0.03 


1 67 +(uw 

2 nn +003 

zuu -0.03 

1 81 +0.05 
1 - 81 -0.04 


1 Q7+0.15 
1 <;7+0.03 

'•-"-o.os 

1 7<r+0.06 
1 -'- , -0.06 

1 77+O.O7 

'•"-o.oe 

7 70+0.78 
'■ ZS -0.62 


0.77/140 
1.04/140 
1.05/628 
1 .Z0/0Z0 

1.17/390 
1.76/390 
0.91/331 
1.09/331 
0.88/330 


8.27 
7.50 
18.40 

10. /U 

9.10 
8.20 
23.30 
21.30 
22.60 


25« 


2001 February 

2004 February 
2004 August 

2005 February 


pow 
diskbb 

pow 
disbb 

pow 
diskbb 

pow 
diskbb 


2 36 +A3S 
Z - JO -0.31 

1 61 + 002 

? 45+O.IO 
Z '^ J -0.10 
1 79+0.06 

'•^-o.oe 

, 6fi +0.12 
1 - DO -0.12 
1 7O+0.08 
1 - z °-0.09 
i 07+0.32 
l -°'-0.29 
1 27+0- 17 


i 5 fi +0.17 


2 11+0.27 
Zil '-0.22 

7 cc+0.09 
z - JJ -0.08 

3-4031 

7 71+0.40 
z - ' '-0.31 


1.14/73 
1.05/73 
1.36/670 
1.18/670 
1.08/310 
0.96/310 
0.97/80 
0.94/80 


7.10 
6.50 
33.30 
31.70 
8.90 
8.50 
8.40 
8.00 


36 


2004 February 
2004 August 


pow+ga 
bb+ga 

pow+ga 
bb+ga 


A fif ;+0.46 
U - DO -0.06 
1 1+ - 20 
a 43+O.O5 

45+0- 10 




Q0+ 008 
u - w -0.08 

0.19^,0.94^04 
0.88!°,$ 

A 1 5+0.01 A A+0.87 

U ' 1;i -0.02' U - u -0.0 


0.86/30 
0.94/30 
0.97/116 
1.23/116 


0.15 
0.16 
0.86 
0.90 


44 


2001 February 

2004 February 
2004 August 

2005 February 


diskbb 
diskbb 
diskbb 
diskbb 


a 21+0.12 
U - ZJ -0.11 
a 21+010 
u ' zl -0.09 
Q21+O.05 
U ' Z1 -0.05 
21+0.18 
u ' zl -0.15 




1 4Q+0.30 

1 30+ 022 
1 - JU -0.17 
1 44+O.IO 

7 it+1.04 
z - -0.52 


0.49/32 
0.73/31 
0.97/60 
1.00/16 


1.40 
0.63 
0.85 
1.50 


52 


2001 February 
2004 February 
2004 August 


pow 
pow 
pow 
RS 


A <j<:+0.18 

73+0-20 
U -' J -0.18 

A CA+0.34 

u ' J? -0.23 


z - s4 -0.79 
3 36 +0 - 63 
Q <;7+0.48 
J ' JZ -0.42 


1 05+ 036 
1 .UJ_ .22 


1.06/15 
1.24/33 
1.18/43 
1.08/42 


0.21 
0.73 
0.30 
0.27 



Spectral variables* 



23 


2001 February 

2004 February 
2004 August 

2005 February 


pow 
pow 
pow 
pow 


o.oo+S 

0.64+0;^ 
0.6££Jg 


a 4 7 +1.72 

1 53 +1 07 
1J -0.46 
1 81+0.58 
1-81-0.49 

2.52! 2 : « 




0.21/8 
0.68/14 
0.79/35 
0.78/15 


0.30 
0.12 
0.20 
0.13 


24 


2001 February 
2004 February 


pow 
pow 


0.88t 22 
0.23i8;}§ 


3 04+0.94 
J - y -0.70 




0.94/38 
1.52/36 


0.09 
0.12 




Fig. 1. — Digitized Sky Survey (DSS) blue band image of (a) the field of view (35 ; x 35') of 
XMM-Newton and (b) the central 10' x 10' of IC342, with detected Chandra X-ray sources over- 
laid. SSS/QSS (diamonds), foreground stars (squares), SNRs (boxcircles) are identified while the 
unclassified X-ray point sources (circles) are also shown. The radii of the unclassified X-ray source 
circles are 15". 
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Table 8 — Continued 



Source ID Obs year XSPEC model N H (10 22 crrT 2 ) r/log(n e O a kT b (keV) X 2 /dof I\ (10 38 erg s" 1 ) 



24 (con't) 




RS 


o.oio;°;8 3 




o.83^;8 7 


1.28/37 


0.08 




2004 August 


pow 


5S+0- 18 
u ""-0.14 


°- JH -0.99 




1.89/53 


0.09 






BB 


1 l+o- 06 

U ' 11 -0.03 




o.i7-8:8i 


1.67/53 


0.09 




2005 February 


pow 


a -i 7+0. 16 
U ' JZ -0.12 






0.95/32 


0.21 


38 


2001 February 


pow+BB+ga 


() c 7 +0.33 


3 07+0.05 
JA "-0.70 


Q<i+0.95 a fi4 +0.23 
u - y;, -0.95> ucn -0.42 


1.12/99 


1.56 






pow+mekal 


n 77+0.17 


2 63+ 018 
Z - OJ -0.16 


24+°' 11 
u - z ^-0.05 


1.09/101 


1.61 




2004 February 


pow+BB+ga 


5<s+0-«> 


7 <ro+0.49 
/ -' ,S -0.63 


n 1 7+0.06 a 77+O.07 

"^'-o.os'^-^-o^s 


1.41/88 


1.32 






pow+mekal 


n 43+0.13 


7 iq+0.25 
z - ly -0.25 


4Q+0.13 


1.15/90 


0.91 




2004 August 


pow+BB+ga 


30+O- 11 
u - 5u -0.08 


1 66+ 034 
i.oo_ 045 


n 9q+0.04 a C9+0.03 
U - Zi -0.04' U ' SZ -0.05 


1.26/113 


1.15 






pow+mekal 


f) ^+0.04 
u ""-0.04 


1 -' u -0.21 


n1 +0 04 


1.18/115 


1.11 




2005 February 


pow+BB+ga 


27+ ' 15 
u - z '-0.15 




Q 7O+0.10 7O+0.49 

u - zo -0.07' u - /y -0.36 


1.26/44 


1.57 






pow+mekal 


A <3Q+0.60 

UjS -0.08 


1.50_j;}f 


<jS+ 008 
UJJ -0.30 


1.28/46 


1.15 



Other bright X-ray sources 



21 


2001 February 

2004 February 
2004 August 

2005 February 


pow 
pow 
pow 

diskbb 
pow 

diskbb 


ojotl 

0-83^ 

o.estli 

34-* 15 
u -^-0.13 
a 7 A+0.84 
u ' /y -0.60 

u - Jy -0.20 


1 77+O.63 
i - // -0.54 

2.00i*; 42 . 

7 77+1.13 
z ' z '-0.87 


7 01 +0.38 
z - -0.29 

1 n+0.96 
1 - 1/ -0.43 


0.81/11 
0.60/14 
0.87/26 
0.62/26 
0.49/12 
0.47/12 


0.48 
0.35 
0.54 
0.50 
0.45 
0.37 


27 


2004 February 
2004 August 


NEI 
NEI 


1.04^ 

o.8i^;H 


y --> -o.oi 

A 1 1+0.01 

v - 11 -0.01 


7 f,4+1.12 

7 1 , +0.56 
z - 11 -0.56 


1.29/14 
0.30/15 


0.07 
0.10 



Note. — This table lists spectral parameters to energy spectra of sources with enough counts (over 100) in a particular observation for fitting. 
The sources in a particular group (i.e. flux, spectral variables, and other bright sources) are listed in ascending order of RA. For the XSPEC 
model, pow refers to power-law; BB refers to blackbody; diskbb refers to disk blackbody; RS refers to Raymond-Smith; ga refers to Gaussian 
lines. The quoted errors are at 90% confidence level as generated by XSPEC. 

"This column shows the photon index (r) and the ionization timescale (n e t) for power-law and NEI models respectively. 
*The temperature refers to blackbody temperature for blackbody, temperature at inner disk radius for disk blackbody, plasma temperature for 
Raymond-Smith and NEI, and line energy in keV for Gaussian. In the case of source 36 and 38, the 2 temperatures are blackbody temperature 
and gaussian line energy in this order. 

c The quoted luminosities are absorbed value. The assumed distance to IC342 is 1.8 Mpc. 
f Sources with significant flux variability (i.e. Sfl ux > 3) 

* Sources with significant spectral variability that exhibit changes in the best-fitted spectral parameters at 90% confidence level. 
4 These two sources are also classified as spectral variables. 




Fig. 2. — Color-color diagram for all sources with more than 20 counts (circles: detected sources; 
squares: spectral models). Also plotted are the theoretical hardness ratios estimated from different 
spectral models. Top to bottom: Power-law model with T = 1.2, 1.7, 2, and 3; Raymond-Smith 
model with £7rs =0.5 keV; Blackbody model with £7bb = 0.1 keV . For each model, varies 
from left as 3 x 10 21 cm" 2 , 5 x 10 21 cm" 2 , and 10 22 cm~ 2 lKongl(l2003k 
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Fig. 3. — The flux variability factor 5fl ux versus the average angle offset from the galactic center of 
IC342. Blue: sources detected in only one observation; Green: sources detected in two observa- 
tions; Magenta: sources detected in three observations; Red: sources detected in all 4 observations; 
Sources with S flux > 3 (above the dash line) are defined to be flux variable. The arrows attached to 
the data point represents the 3 a lower limit of the variability factor. 
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Fig. 4. — Long-term light curve of the X-ray sources in IC342: source 17 (top-right), source 19 
(top-left), source 25 (bottom-right), and source 38 (bottom-right). Filled circle represent Chandra 
HRC-I observations, triangles represent XMM-Newton EPIC observations, and diamonds represent 
ROSAT HRI observations. Errors for XMM-Newton data points are la gaussian limits from XSPEC 
while for Chandra and ROSAT data points are derived poisson errors of the source counts. It 
was noted that the XMM-Newton and ROSAT fluxes of source 38 could suffer from background 
contamination (refer to Section l5.2.3l) . All luminosities are unabsorbed values scaled to a distance 
of 1.8 Mpc. 




Fig. 5. — Hardness-intensity diagram for all sources detected in the four XMM-Newton observa- 
tions. All HR1, HR2, and luminosity are average values as given in Tabled 
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Fig. 6. — XMM-Newton EPIC-PN energy spectra of the three brightest X-ray sources: source 19 
(IC342 X-l; top-left), source 25 (IC342 X-2; top-right), and source 38 (IC342 X-3 nuclear source; 
bottom-right). The fourth plot is the source 27 (SNR candidate; bottom-left). The 2001 February 
(black triangle), 2004 February (red cross), 2004 August (green circle), and 2005 February (blue 
star) data are shown together with the lines showing the best-fit spectra. The best-fit parameters 
are given in Table [8] 



